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Introduction 


X-ray  mammography  is  currently  the  clinically  accepted  modality  for  breast 
eaneer  sereening.  Mammography,  however,  is  particularly  unreliable  and  recent  trials 
have  shown  that  70-80%  of  mammographically  indeterminate  lesions  that  progress  to 
surgical  biopsy  are  benign.  Development  of  an  accurate,  non-invasive  imaging  test 
with  low  false  negative  rates  would  allow  reduction  of  the  two  step  process  of  surgieal 
biopsy  followed  by  surgical  lesion  removal  to  a  single  step,  as  well  as  providing  a  tool  for 
guiding  treatment  deeisions.  It  may  be  possible  to  achieve  this  through  the  development 
of  a  non-invasive  technique  that  measures  microvessel  density.  Tumour  progression  is 
known  to  be  dependent  on  the  ability  to  stimulate  the  growth  of  new  blood  vessels  that 
supply  nutrients  and  oxygen  to  the  expanding  tumour.  Based  on  recent  studies  that 
indicate  that  the  onset  of  angiogenesis  ean  oceur  before  a  tumour  reaehes  its  maximum 
diffusion-limited  growth  potential(l-2  mm  diameter)  and  may  oeeur  well  before  the 
switch  to  a  malignant  phenotype,^"’^  it  has  been  hypothesized  that  the  eharaeteristie  new 
vessel  growth  associated  with  tumour  angiogenesis  should  allow  differentiation  between 
malignant  and  benign  breast  lesions. The  purpose  of  this  research  is  to  develop  a 
non-invasive  predictor  of  malignancy  through  the  design  and  preliminary  validation  of  a 
new  form  of  magnetic  resonance  imaging  (MRI)  called  the  IntraVoxel  Ineoherent  Motion 
(IVIM)  method  which  utilises  high  main  magnetie  field  strength  and  ultra-strong 
magnetic  field  gradients.  This  teehnology  will  permit  completely  non-invasive, 
quantitative  charaeterization  of  tumour  mierovasculature  at  high  spatial  resolution  over 
the  entire  breast. 
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Body 

Hardware  Development 

A  high  strength  gradient  coil  has  been  designed.  The  prototype  coil  was  built 
with  an  inner  diameter  of  5cm,  which  is  ideal  for  mice.  A  second  coil  has  now  been  built 
that  is  more  suitable  for  rats.  This  coil  has  an  inner  diameter  of  24cm,  which  can 
accommodate  a  rat  plus  the  equipment  used  to  monitor  the  animal.  The  maximum 
gradient  strength  achievable  with  this  coil  is  300mT/m,  which  is  almost  ten  times  greater 
than  typical  clinical  systems  and  is  sufficient  for  the  high-resolution  imaging  experiments 
proposed  in  this  grant.  Testing  of  the  coil  on  both  a  l.STesla  (T)  clinical  MRI  system  and 
a  4T  research  MRI  system  is  under  way.  It  was  found  that  wire  density  and  power 
dissipation  considerations  required  that  the  coil  be  constructed  in  a  multilayer 
configuration.  Details  regarding  the  design  and  construction  of  the  coil  are  given  in  the 
paper  in  Appendix  I.  This  paper  has  been  accepted  by  the  journal  MAGMA. 

Pulse  Sequence  Development 

The  IVIM  pulse  sequence  has  been  designed  based  on  a  standard  diffusion- 
weighted  imaging  sequence.  The  sequence  consists  of  a  spin  echo,  echo  planar  imaging 
sequence  with  the  addition  of  strong  magnetic  field  gradients  that  cause  the  contrast  in 
the  images  to  be  weighted  based  on  the  diffusion  of  water  in  the  microvascular  and  extra- 
vascular  compartments.  We  have  incorporated  a  simple  phase  correction  method  that 
allows  post-acquisition  correction  of  the  images  for  small  scale  motion  of  the  object 
being  imaged. 

A  computer  program  has  been  'written  in  an  interactive  windows  environment  that 
allows  analysis  of  the  diffusion  weighted  images  obtained  with  the  IVIM  pulse  sequence. 
The  program  allows  calculation  of  the  diffusion  coefficients  either  for  each  pixel  in  the 
images,  or  for  the  average  values  in  a  region  of  interest.  Validation  of  the  program  has 
been  performed  using  simulated  data.  Further  validation  of  both  the  program  and  the 
IVIM  sequence  has  been  performed  using  phantoms  that  are  spherical  containers  filled 
-with  water  and  acetone.  The  results  obtained  (see  Appendix  II)  agree  well  with  those 
found  in  the  literature. 
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Validation  of Microvascular  Imasins 

A  tumour  model  has  been  implemented  in  rats  using  N-ethyl-N-nitrosourea 
(ENU).  ENU  is  known  to  be  a  potent  carcinogen  that  induces  well  vascularized 
mammary  tumours  in  rats?^  This  tumour  model  was  chosen  because  a  range  of  benign 
and  malignant  tumours  develop  depending  on  the  dose  of  ENU  administered  and  because 
the  tumour  morphology  and  angiogenic  behavior  closely  mimic  breast  tumours  in 
humans.  Tumours  have  been  successfully  grown  in  30  rats  that  were  given  a  high  dose 
injection  of  ENU  and  in  6  rats  given  a  low  dose.  IVIM  images  of  the  tumours  have  been 
obtained  for  each  rat.  High  signal  to  noise  ratios  (SNR)  are  required  in  the  images  in 
order  to  observe  the  microvascular  component  of  the  measured  signal.  SNR’s  in  the 
images  have  been  at  least  400,  which  should  be  sufficient  for  this  study.  Two  of  the  rats 
were  sacrificed  immediately  after  imaging,  and  imaging  was  repeated  on  the  sacrificed 
animal  in  order  to  verify  that  the  IVIM  component  disappears  after  cessation  of  blood 
flow  and  is  in  fact  due  to  microvascular  blood  flow.  Preliminary  results  (given  in 
Appendix  III)  show  a  fast  decay  component.  The  magnitude  of  the  diffusion  coefficient 
for  this  component  is  in  on  the  order  of  that  expected  for  the  microvascular  blood  flow. 
This  fast  decay  component  disappears  after  cessation  of  blood  flow.  These  results 
suggest  that  the  IVIM  method  is  capable  of  measuring  the  tumour  microvasculature. 
Further  analysis  of  the  acquired  images  is  underway.  The  tumours  have  been  resected 
and  pathology  will  be  performed  on  each  tumour  to  determine  the  degree  of  malignancy 
and  the  density  of  microvessels  in  the  tumours.  The  pathological  results  will  be 
correlated  with  the  MR  images. 
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Key  Research  Accomplishments 


•  Design  of  a  mouse-sized  gradient  coil  with  gradient  strengths  that  exceed  2000mT/m. 
(Approximately  50  to  100  times  higher  than  typical  clinical  MRI  systems.) 

•  Design  of  a  gradient  coil  with  an  inner  diameter  of  24cm,  which  is  large  enough  to 
accommodate  a  rat  and  animal  monitoring  equipment.  Maximum  gradient  strengths 
achieved  with  this  coil  will  be  300mT/m. 

•  Implementation  of  a  fast  MR  pulse  sequence  that  allows  acquisition  of  diffusion 
weighted  images  with  high  signal  to  noise. 

•  Design  and  coding  of  a  computer  program  that  will  allow  automated  analysis  of  the 
diffusion  weighted  images  for  the  calculation  of  diffusion  coefficients. 

•  Validation  of  the  MR  pulse  sequence  analysis  program  using  simulated  data  and 
images  acquired  from  substances  with  known  diffusion  coefficients. 

•  Implementation  of  a  breast  tumour  model  in  rats  that  produces  tumours  with  a  range 
of  malignancies  from  benign  to  malignant. 

•  Acquisition  of  in  vivo  diffusion  weighted  images  of  the  tumours  in  the  rats  using  the 
IVIM  method,  and  preliminary  analysis  of  the  images. 
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Reportable  Outcomes 


Presentations  and  Publications 

1 .  The  design  and  construction  of  the  gradient  coils  resulted  in  two  peer-reviewed 
publications^^’^'*  and  two  conference  presentations.^^’^^  These  two  publications  and  a 
third  publication,  recently  accepted  by  MAGMA,  are  appended  (Appendices  I,  IV, 
andV). 

2.  The  preliminary  results  of  the  acquisition  and  analysis  of  the  images  acquired  from 
the  animal  tumour  model  were  presented  at  the  Era  of  Hope  Meeting,  Atlanta,  Jime  8- 
12, 2000. 

Funding 

Two  post-doctoral  fellowships  have  been  awarded  to  perform  the  research 
supported  by  this  grant: 

1 .  To  Dr.  Paula  Gareau  by  the  Department  of  Defense  Breast  Cancer  Research  Program 
of  the  United  States  Army  Medical  Research  and  Materiel  Command.  Grant  Number 
BC980729. 

2.  To  Dr.  Lanette  Friesen-Waldner  by  the  Natural  Sciences  and  Engineering  Research 
Council  of  Canada. 

Dr.  Lanette  Friesen-Waldner  has  recently  applied  for  a  postdoctoral  traineeship  to 
the  Department  of  Defense  U.S.  Army  Medical  Research  and  Materiel  Command’s  2000 
Breast  Cancer  Research  Program.  The  research  proposed  in  this  application  is  based  on 
the  work  supported  by  this  award. 
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Conclusions 


We  have  designed  a  mouse-sized  gradient  coil  and  a  gradient  coil  large  enough  to 
accomodate  rats  and  animal  monitoring  equipment.  Both  coils  achieve  gradient  strengths 
significantly  higher  than  typical  clinical  MRI  systems.  Based  on  wire  density  and  power 
dissipation  considerations,  we  have  concluded  that  small  coils  with  high  gradient 
strengths,  such  as  those  required  for  high  resolution  IVIM  imaging,  require  multilayer 
configurations.  The  design  for  the  coils  constructed  here  should  be  scalable  to  a  breast¬ 
sized  coil.  Such  a  coil  will  allow  imaging  of  the  human  breast  at  resolutions  that,  until 
now,  have  been  unattainable,  thus  allowing  visualisation  of  very  small  lesions  within  the 
breast. 


An  MR  pulse  sequence  has  been  implemented  on  our  1.5T  clinical  MRI  system 
that  will  allow  us  to  acquire  IVIM  images  with  high  signal  to  noise  ratios  (>400).  A 
computer  program  has  been  written  to  analyse  the  IVIM  images.  The  pulse  sequence  and 
analysis  program  have  been  validated  using  substances  with  known  diffusion 
coefficients.  Images  have  also  been  acquired  of  tumours,  which  were  developed  using  an 
animal  tumour  model  that  we  have  implemented.  Analysis  of  the  images  has  begun  and 
the  preliminary  results  suggest  that  the  tools  developed  here  should  be  able  to  measure 
non-invasively  the  microvasculature  in  breast  tumours.  This  would  allow  diagnosis  of 
breast  lesions  in  vivo  and  would  also  provide  a  means  of  following  the  effects  of 
treatments,  particularly  treatment  with  anti-angiogenic  pharmaceuticals;  however,  further 
testing  of  our  methods  is  required  before  definitive  conclusions  can  be  drawn. 
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Abstract: 


There  is  great  interest  in  the  non-destructive  capabilities  of  magnetic  resonance  microscopy 
for  studying  murine  models  of  both  disease  and  normal  function;  however,  these  studies  place 
extreme  demands  on  the  MR  hardware,  most  notably  the  gradient  field  system.  We  designed, 
using  Constrained  Current  Minimum  Inductance  methods,  and  fabricated  a  complete  three-axis 
gradient  coil  set  that  utilizes  interleaved,  multilayer  axes  to  achieve  maximum  gradient  strengths 
of  over  2000mT/m  in  rise  times  of  less  than  SOjxs  with  an  iimer  coil  diameter  of  5cm.  The  coil 
was  wire- wound  using  a  rectangular  wire  that  minimizes  the  deposited  power  for  a  given  gradient 
efficiency.  Water  cooling  was  also  incorporated  into  the  coil  to  assist  in  thermal  management. 
The  duty  cycle  for  the  most  extreme  cases  of  single  shot  EPI  is  limited  by  the  thermal  response 
and  expressions  for  maximum  rates  of  image  collection  are  given  for  burst  and  continuous  modes 
of  operation.  The  final  coil  is  capable  of  the  collection  of  single  shot  EPI  images  with  6mm  field 
of  view  and  94pm  isotropic  voxels  at  imaging  rates  exceeding  50s‘^ 
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Introduction 


High  speed  magnetic  resonance  microscopy  in  murine  models  of  human  disease  and  fiinction 
is  of  fundamental  importance  in  a  range  of  biomedical  research  areas  ( 1 ,2).  Non-destructive  MR 
histology”  of  brain  lesions  in  small  animals  (3,4),  studies  of  the  biophysical  basis  of  fMRI  in  a 
murine  model  (5),  in-vivo  diffusion-tensor  microscopy  (6),  microscopic  studies  of  in  utero  murine 
embryos  (7),  in-vivo  microscopic  studies  of  the  murine  carotid  artery  (8),  and  murine  cardiac  MRI 
(9,10)  are  all  active  areas  of  research  requiring  in  vivo  magnetic  resonance  microscopy.  Cardiac 
MRI  in  the  mouse  heart  is  die  most  technically  challenging  because  of  the  combination  of  small 
organ  size  (<8mm)  and  very  fast  heart  rate  (>450  beats  per  minute).  The  technical  demands 
placed  on  the  system  hardware  in  these  applications  are  extreme,  particularly  the  gradient  coil 
requirements.  The  small  size  of  the  mouse  is  of  course  an  advantage  in  that  it  allows  the  use  of 
small,  customized  gradient  coils.  This  paper  discusses  the  special  considerations  involved  in 
designing  such  mouse-specific  coils,  and  details  the  design,  construction,  and  imaging  perfor- 
mance  of  one  specific  design. 

AboVe  all  else,  very  high  gradient  strengths  are  required  to  achieve  the  high  spatial  and  tempo- 
ral  resolution  necessary  for  advanced  MRI  of  small  animals.  Using  the  fundamental  relationship 
between  gradient  strength,  sampling  rate,  and  image  field  of  view  (FOV)  (1 1),  it  is  easily  found 
that  a  6mm  FOV  requires  a  gradient  strength  of  2000mT/m  when  using  a  receive  signal  sweep 
width  of  500kHz,  typical  of  the  maximum  sampling  rate  of  recent  generation  MRI  systems.  This 
gradient  strength  is  approximately  two  orders  of  magnitude  greater  than  that  produced  by  the 
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whole-body  gradients  of  current  clinical  scanners;  therefore,  it  is  clear  that  a  customized  gradient 
coil  is  necessary  for  high  quality  mouse  organ  imaging. 

For  any  given  amplifier  and  gradient  coil  combination,  the  maximum  rise  time  of  the  gradient 
waveforms  is  proportional  to  the  coil  inductance  (L).  For  this  reason,  coils  are  often  designed 
with  minimum  inductance  for  a  given  gradient  strength.  This  minimizes  the  total  time  spent 
switching  in  any  gradient  waveform  train,  thereby  minimizing  TR,  TE,  and  other  critical  timing 
parameters.  The  switching  speed  requirements  for  high  field  mouse  imagers  are  especially 
extreme.  For  example,  at  7T  and  above,  the  T2*  of  murine  myocardial  tissue  is  9ms  or  less  and 
the  T2*  of  gray  and  white  matter  are  18ms  and  12ms,  respectively  (12,13).  Acquiring  a  64-lobe 
EPI  train  within  20ms,  using  500kHz  flat-top-only  sampling  and  a  6mm  FOV,  will  require  a  mini¬ 
mum  efficiency  of 2000mT/m  and  a  maximum  inductance  of  200|j.H.  This  inductance  limit  is  up 
to  an  order  of  magnitude  lower  than  that  of  typical  clinical  gradient  coils. 

A  final  and  often  overlooked  consideration  is  the  gradient  coil  duty  cycle,  determined  by  the 
thermal  response  of  the  coil,  which  is  in  turn  largely  determined  by  the  resistance  of  the  coil  load. 
Because  many  researchers  will  wish  to  connect  their  mouse-sized  gradient  coil  to  a  normal  clini¬ 
cal  gradient  driver  system,  the  coil  must  be  capable  of  carrying  high  rms  currents  (in  the  range  of 
lOOA).  For  resistances  between  0.5  and  lO,  this  corresponds  to  power  deposition  of  5  to  lOkW  in 
a  physically  small  coil  (low  thermal  mass)  and  there  is  the  potential  for  a  significant  heating  prob¬ 
lem  and  a  limited  duty  cycle.  Because  of  the  small  coil  size  needed  to  achieve  the  strength  and 
speed,  the  area  available  for  wire  cross-section  is  greatly  reduced.  In  the  following  sections  it  will 
be  demonstrated  that  simply  scaling  down  existing  body  gradient  designs  is  insufficient  when  this 
power  deposition  issue  is  considered.  The  way  to  achieve  the  required  performance  is  to  move  to 
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multilayer  coil  designs  while  using  rectangular  wire  to  maximize  copper  cross-section.  This 
method  in  turn  creates  the  need  to  specifically  consider  an  interleaving  strategy  for  organizing  the 
many  layers  of  the  coil  into  a  complete  three  axis  design.  We  outline  a  design  strategy  for  inter¬ 
leaved  multilayer,  rectangular  wire  gradient  coils  in  this  paper,  and  show  the  expenmental  perfor¬ 
mance  of  a  three-axis  mouse-optimized  coil  created  using  this  strategy. 

Theory 

Gradient  Coil  Scaling 

As  was  noted  in  the  introduction,  the  advantage  of  a  customized  mouse  gradient  coil  insert  is 
that  it  can  be  built  to  a  small  radius  just  large  enough  to  encompass  the  mouse  and  RF  coil  assem¬ 
bly.  We  decided  to  design  our  prototype  coil  to  an  inner  physical  radius  of  2.5cm,  and  therefore 
the  innermost  winding  layer  of  the  coil  would  have  a  radius  of  approximately  3cm.  This  is  one 
order  of  magnitude  smaller  than  would  be  typical  for  normal  whole  body  gradient  coils.  The  first 
question  we  sought  to^answer  was  what  sort  of  performance  could  be  expected  simply  by  scalmg 
a  human  whole-body  gradient  coil  down  to  mouse  size. 

The'  four  most  important  quantities  defining  the  performance  of  a  gradient  coil  are  the  effi¬ 
ciency  (ti),  inductance  (L),  resistance  (R),  and  temperature  rise  (AT)  for  a  given  input  current. 
The  general  laws  for  scaling  these  quantities  from  values  at  radius  ao  and  peak  winding  density 
(defined  as  the  number  of  wires  per  cm  used  to  wind  the  coil  at  the  densest  region  of  the  wire 

pattern)  to  new  radius  a  and  density  are. 


T1  =  fio- 


[1] 
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[2] 


R  =  Ro- 


[3] 


where  is  the  cross-sectional  area  of  the  wire.  Resistance  is  of  importance  in  this  application 
because  of  its  influence  on  the  thermal  response  of  the  coil.  The  average  temperature  increase  of 
a  coil  is  proportional  to  resistance  and  inversely  proportional  to  the  total  volume  of  the  coil  wall 
(assuming  that  the  material  properties  remain  constant  with  temperature).  Assuming  that  the  vol¬ 
ume  scales  as  the  cube  of  the  radius  and  that  A^  scales  inversely  with  the  square  of  the  winding 
density  (as  they  would  for  uniformly  scaled  coil  and  wire  dimensions),  then  the  temperature 
change  at  any  constant  rms  current  is  expected  to  scale  as: 


AT  =  ATq  • 


[4] 


'V 

Figure  1  contains  plots  of  r\,  L,  and  AT  as  a  function  of  radius  for  several  fixed  values  of  winding 
density.  The  starting  point  used  to  generate  these  plots  was  a  hypothetical  “typical”  body-sized 
gradient  coil  design  with  a  radius  of  30cm,  efiiciency  of  0.12mT/ni/A,  inductance  of  1200|j,H, 
resistance  of  SOOmH,  temperature  change  of  5°C  at  an  rms  current  of  lOOA,  and  a  maximum 
winding  density  of  2cm*  ^ . 


For  a  clinical  amplifier  system  capable  of  delivering  a  peak  current  of  200 A,  the  gradient  effi¬ 
ciency  required  to  achieve  2000mT/m  is  lOmT/m/A.  Examining  figure  la  we  note  that  a  winding 
density  of  16cm''  or  higher  is  required  to  achieve  lOmT/m/A  for  a  radius  of  3cm.  This  is  itself  a 
limitation  in  that  constructing  such  a  dense  wire  pattern  is  quite  difficult  with  standard  methods; 
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furthermore,  the  wire  would  have  to  be  of  such  small  cross  section  as  to  preclude  the  high  peak 
currents  (on  the  order  of  the  200A)  required.  Construction  methods  developed  in  our  lab  cur¬ 
rently  produce  a  maximum  wire  density  of  approximately  8cm’ ‘  on  a  single  transverse  layer,  but 
according  to  the  efficiency  plot  this  only  produces  approximately  half  of  the  required  strength. 
One  way  to  reach  higher  effective  density  wire  patterns  with  lower  density  fabncation  methods  is 
to  divide  the  winding  of  a  single  layer  across  two  or  more  layers,  each  individual  layer  wound  at 
the  highest  density  achievable. 

The  inductance  variation  depicted  in  figure  lb  indicates  that  coils  at  the  highest  values 
considered  fall  below  the  200pH  upper  limit  allowed  for  this  application.  The  figure  mdicates  that 
a  single  layer  design  wound  at  8cm’‘  would  have  an  inductance  of  20pH  at  a  3cm  radius.  To  a 
first  approximation,  the  inductance  of  a  double  layer  8cm-‘  coil  operated  in  series  is  equivalent  to 
that  of  an  equal  radius  single  layer  design  wound  to  16cm’\  i.e.  80p,H  at  a  3cm  radius;  therefore, 
it  is  expected  that  the  two  layer  configurations  needed  to  achieve  the  efficiencies  discussed  m  the 
previous  paragraph  wj|l  also  result  in  inductance  values  below  the  specified  upper  limit. 

The;  thermal  response  of  small,  high  winding  density  coils  represents  the  dominant  practical 
limitation,  as  is  evident  from  figure  Ic.  Scaled  to  the  desired  mouse  coil  radius,  all  winding  den¬ 
sities  result  in  steady-state  coil  temperatures  orders  of  magnitude  higher  than  would  be  tolerable. 
For  this  reason,  we  expect  that  all  mouse  coils  meeting  the  specifications  outlined  will  have 
extreme  duty  cycle  limitations  unless  special  effort  is  made  to  minimize  the  temperature 
increases.  We  argue  in  this  paper  that  this  consideration  leads  to  a  dominant  requirement  for  min¬ 
imum  resistance  configurations,  which  requires  the  combination  of  rectangular  cross-section  wire 
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with  multilayer  designs  in  order  to  achieve  the  smallest  resistance  for  a  given  effective  winding 
density. 


Heating 

The  temperature  response  of  a  gradient  coil  when  driven  by  an  rms  current  (Inns)  can  be 
approximated  by  an  exponential  approach  to  an  equilibrium  temperature  (14)  and  the  following 
equation  relating  temperature  (T)  and  time  (t)  applies  to  both  heating  and  cooling  of  the  coil; 


T  = 


To  +  at, 3 


[5] 


where  ATgg  is  the  steady  state  temperature  change  at  the  measurement  location,  T^  is  the  tempera¬ 
ture  at  the  start  of  the  measurement,  and  x  is  a  time  constant  characteristic  of  the  axis  and  the 
degree  of  cooling  applied.  For  the  heating  phase,  a  steady-state  heating  parameter  can  be  defined; 


K  = 


AT 


SS 


iLr 


16] 


where  is  the  rms  current  applied  and  R  is  the  resistance  of  the  axis.  To  a  first  approximation, 
ATss  is  simply  proportional  to  the  input  power  while  x  is  independent  of  it  (14);  therefore,  assum¬ 
ing  that  resistance  is  constant  over  the  range  of  temperatures  in  the  experiment,  k|,  is  a  constant  for 
each  axis  and  rate  of  water  cooling  and  it  can  be  used  to  calculate  the  steady-state  temperature 
increase  for  any  input  power. 


The  fact  that  x  is  independent  of  the  input  power  means  that  it  is  an  operationally  independent 
parameter  that  can  be  used  to  define  the  transition  between  two  modes  of  coil  operation.  These 
two  modes  require  individual  consideration  when  calculating  duty  cycles.  To  achieve  scan  repeti- 
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tion  (TR)  times  much  less  than  x,  the  coil  may  be  operated  in  a  burst  mode,  during  which  it  is  run 
at  high  power  for  a  short  time  (tj,,  the  burst  time)  and  allowed  to  heat-up  at  an  unsustainable  rate 
(figure  2a).  This  would  then  have  to  be  followed  by  a  dead  time,  t^,  during  which  the  coil  is 
allowed  to  cool  sufficiently  before  starting  another  burst  phase.  This  operation  mode  would  be 
required  when  imaging  at  very  high  rates  is  important  but  required  only  for  short  blocks  of  time; 
an  example  might  be  cardiac  imaging.  For  this  mode,  the  maximum  number  of  single  shot  EPI 
images  that  can  be  collected  during  each  burst  phase  was  calculated  as  a  function  of  the  pulse 
sequence  and  coil  thermal  parameters  such  that  the  coil  temperature  increases  were  maintained  at 

tolerable  levels. 

For  scan  times  on  the  order  of  x  and  longer,  the  coil  must  be  run  in  a  sustainable  fashion,  and 
the  steady-state  thermal  response  is  of  primary  importance.  In  this  continuous  mode  (figure  2b), 
the  minimum  TR  for  repeated  single  shot  EPI  imaging  was  calculated  as  a  function  of  the  pulse 
sequence  parameters  such  that  the  corresponding  rms  current  resulted  in  a  tolerable  steady-state 
temperature  increase.  Burst  and  continuous  mode  operation  are  analyzed  in  more  detail  in  the  fol¬ 
lowing  paragraphs. 

The  purpose  of  burst  mode  operation  is  to  allow  images  to  be  collected  at  the  minimum 
sequence  repetition  time  defined  by  timing  rather  than  by  heating  considerations.  The  length  of 
time  that  this  can  be  sustained  (tb)  is  determined  by  the  rms  current  applied  during  the  burst 
(In^  b)»  *e  temperature  excursion  allowed  to  develop  during  the  burst  (ATfc),  and  the  heating 
parameters  of  the  coil.  Ignoring  the  heating  contributions  from  the  phase  encode  and  slice  select 
axes,  the  rms  current  that  exists  during  each  single  shot  EPI  collection  can  be  approximated  as; 
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1/2 


[7] 


^j^27c-SW^  K  (2  Lr  27t-SW^  ^v] 
rms_b  -  l^yti^FO V  J  ■  VtR  ■  [3  ■  V  ■  y ti,FO V,  S W  J 

where:  SW  is  the  receive  signal  sweep  width,  y  is  the  gyromagnetic  ratio  of  the  nucleus,  r\j.  is  the 
gradient  efficiency  of  the  readout  axis,  FOVj.  is  the  field  of  view  in  the  direction  of  the  readout 
axis,  Np  is  the  number  of  phase  encode  lines,  TR  is  the  time  between  EPI  collections,  is  the 
inductance  of  the  readout  axis,  V  is  the  voltage  applied  during  gradient  slewing,  and  Nj.  is  the 
number  of  readout  points  collected.  Applied  for  an  infinite  amount  of  time,  this  rms  current 
would  typically  result  in  an  extremely  high  steady  state  temperature  increase  found  by  solving 
equation  [6]  for  the  temperature: 

AT,,_b  =  k|,.lL.b-R'  [8] 

This  steady  state  temperature  increase  is  used  within  equation  [5]  to  obtain  the  coil  temperature  as 
a  fimction  of  time.  For  burst  mode  operation,  the  initial  slope  of  the  temperature  response  is  criti¬ 
cal  and  the  temperature  cxirve  can  be  approximated  for  tj, «  x  to  give: 


where  the  temperature  increase  at  any  given  instant  is  defined  as  ATj,  =  T  -  Tq.  The  maximum 
number  of  images  (Nburst)  can  be  collected  within  a  single  burst  is  simply  the  burst  time 
divided  by  the  minimum  sequence  TR.  Solving  for  tj,  via  equations  [9]  and  [10]  and  dividing  by 
TR  we  obtain: 
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Equations  [7]  and  [10]  yield  the  following  general  expression  for  the  maximum  number  of  images 
that  may  be  collected  as  a  function  of  the  pulse  sequence  parameters,  the  maximum  tolerable  ATb, 
and  the  readout  axis  heating  parameters  (kh  and  x): 


ATb  X  ^yilrFQVrf  .R  . b .  ‘  [11] 

Nburst  Np  'kbR'^Ti-SW  J  [3  V  yilrFOV,  SWj 

Using  this  expression,  Nbux.t  was  computed  as  a  function  of  the  FOV  for  different  imaging  param¬ 
eters  and  choice  of  gradient  axis  used  for  readout.  Note  that  TR  has  factored  out  of  the  expres¬ 
sion.  Although  shorter  TRs  allow  for  faster  imaging,  they  also  result  in  shorter  allowable  tb*, 
therefore,  the  two  effects  cancel  each  other  and  the  number  of  images  allowed  stays  constant.  The 
last  piece  of  data  that  is  necessary  to  characterize  burst  mode  operation  is  the  length  of  dead  time 
(td)  required  between  bursts.  One  simple  specification  is  to  require  that  the  rms  current  calculated 
over  an  entire  cycle  +  td)  lead  to  an  rms  temperature  increase  (AT^ns)  over  baseline  equal  to 
ATb,  the  burst  time  temperature  increase  (see  Figure  2a).  In  this  case,  the  following  simple  rela- 
tionship  between  and  x  can  be  shown  to  hold. 


Equation  [13]  indicates  that  imaging  can  only  be  done  for  a  period  of  tfe  once  every  heating  time 
constant  of  the  readout  axis.  The  long-term  average  rate  of  image  collection  (i.e.  averaged  over 
many  burst  cycles)  is  independent  of  x  because  although  a  smaller  x  allows  fewer  images  to  be 
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collected  per  burst  (equation  [1 1]),  it  also  allows  more  burst  periods  per  unit  time  (equation  [13]) 
and  the  two  effects  cancel. 


Continuous  mode  operation  is  simply  governed  by  adjusting  the  single  shot  EPI  image  repeti¬ 
tion  time,  TR,  such  that  the  coil  can  be  sustainably  operated.  A  tolerable  steady  state  temperature 
excursion,  AT^,  must  be  chosen  (figure  2b),  and  equation  [6]  is  used  to  calculate  the  correspond¬ 
ing  Iniis  current.  Equation  [7]  is  solved  for  TR  and  combined  with  equation  [6]  to  yield: 


=  N. 


'mm 


f  k,.'. 

I  pn-SWV 

‘4- 

Lr 

UtJ 

I  lyilrFOvJ 
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V 

271  •  SW  1 
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[13] 


Using  tiiis  expression,  minimum  TR  was  computed  as  a  function  of  the  FOV  for  different  imaging 
parameters  and  choice  of  gradient  axis  used  for  readout. 


Methods 

Optimized  Layer  Design 

With  the  need  for  multilayer,  rectangular  wire  designs  identified,  there  remains  the  issue  of 
finding'the  best  design  for  the  individual  layers.  The  three  axes  of  a  gradient  coil  set  can  be  bro¬ 
ken  into  two  design  groups:  the  transverse  (X  and  Y)  and  the  longitudinal  (Z)  axes.  The  X  and  Y 
axes  use  the  same  design,  one  obtained  from  the  other  via  a  ninety  degree  rotation  about  the  longi¬ 
tudinal  axis.  The  strategy  for  selecting  the  optimal  design  is  the  same  for  both  transverse  and  lon¬ 
gitudinal  coils  and  we  will  focus  on  the  former. 

All  designs  were  obtained  using  the  method  of  Constrained  Current  Minimum  Inductance 
(CCMI)  reported  previously  (15).  Specifically,  the  CCMI  method  is  a  solution  for  the  minimum 
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inductance  current  density  satisfying  a  set  of  field  constraints,  and  subject  to  a  set  of  constraints 
on  the  current  density  itself.  The  current  constraints  define  the  region  over  which  the  current  is 
allowed  to  flow  and  therefore  specify  the  design’s  aspect  ratio,  defined  as  the  ratio  of  the  coil  s 
total  length  to  its  diameter.  The  field  constraints  are  placed  in  such  a  way  as  to  define  the  region 
of  uniformity  (ROU).  By  varying  both  the  arrangement  of  the  field  constraints  and  the  lunits  of 
the  current  constraints,  designs  can  be  obtained  for  any  combination  of  ROU  and  aspect  ratio,  and 
for  each  design  the  gradient  efficiency,  inductance  and  winding  density  can  be  calculated.  The 
simplicity  and  computational  efficiency  of  the  CCMI  method  offer  the  unique  opportunity  to 
directly  model  all  possible  designs  in  the  ROU,  aspect  ratio  parameter  space.  By  plotting  design 
performance  figures  over  this  space  at  the  winding  density  and  radius  of  the  desired  coil,  the  entire 
range  of  achievable  performance  is  made  apparent  while  trade-offs  between  imagmg  region  size 
and  coil  size  can  be  rapidly  visualized  and  understood.  The  choice  of  the  optimal  design  then 
becomes  a  straightforward  matter  of  locating  the  best  coordinate  in  this  space  for  the  given  appli¬ 
cation. 

An  additional  complexity  in  our  application  was  the  requirement  for  a  small  side  aperture 
(1cm  iri  diameter)  for  lateral  sample/probe  access.  This  feature  required  the  application  of  current 
constraints  (via  CCMI)  about  the  z  centre  of  the  design  and  the  effects  of  this  process  will  be 
described  in  more  detail  elsewhere. 

Interleaving 

In  the  theory  section  it  was  noted  that  a  multilayer  configuration  would  be  needed  in  this 
application  to  achieve  the  required  gradient  strengths.  It  was  also  suggested  that  rectangular  wire 
would  be  of  significant  benefit  in  minimizing  the  resistance  for  any  winding  density,  and  therefore 
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maintaining  tolerable  thermal  response  for  the  eoil.  The  complicating  issue  in  small  radius,  mul¬ 
tilayer  coils  with  rectangular  wire  is  that  the  thickness  of  each  winding  layer  becomes  a  signifi¬ 
cant  fraction  of  the  radius.  This  is  not  a  major  concern  when  only  one  axis  of  a  gradient  coil  is 
being  built;  however,  it  becomes  a  problem  when  all  axes  of  the  coil  are  being  considered 
together.  In  most  imaging  applications  it  is  not  often  useful  to  have  one  axis  significantly  stronger 
than  the  others.  If  a  multilayer  coil  is  simply  wound  consecutively  from  the  first  axis  to  the  last  as 
in  figure  3a,  the  width  of  the  layers  will  force  the  effective  radius  of  the  outermost  axis  to  be  sig¬ 
nificantly  larger  than  the  innermost,  resulting  in  an  appreciable  performance  imbalance.  The 
alternate  approach  is  to  interleave  the  layering  order  in  such  a  way  as  to  equalize  the  strength,  an 
example  of  which  is  shown  in  figure  3b.  The  interleaving  question  becomes  that  of  determining 
the  effects  of  different  distributions  of  the  layers  of  the  three  axes  into  the  allowed  radial  posi¬ 
tions. 

In  preliminary  fabrication  tests  for  this  design,  it  was  determined  that  an  entire  single  layer, 
including  return  paths'  and  potting  considerations,  could  be  achieved  with  a  wall  thickness  of 
7mm  With  the  radius  of  the  inner  bore  set  at  2.5cm  and  assuming  a  five  layer  coristruction  (the 
choice  of  a  five  layer  construction  is  explained  in  the  results  section),  the  outer  radius  of  the  coil 
would  be  6.0cm.  The  effective  wire  pattern  radius  for  each  layer  was  chosen  to  be  2.75mm  larger 
than  the  irmer  radius  of  that  layer,  yielding  radii  of  2.775cm,  3.475cm,  4.175cm,  4.875cm,  and 
5.575cm  for  layers  one  through  five,  respectively.  There  are  fifteen  distinct  possibilities  for  dis¬ 
tributing  the  three  axes  over  the  five  layers  and  they  were  each  modeled  by  scaling  the  chosen 
transverse  and  longitudinal  coil  designs  to  the  above  radii.  The  winding  density  was  set  to  8cm  ^ 
for  each  of  the  transverse  layers  and  to  lOcm"^  for  the  longitudinal  layer.  The  two  layers  of  each 
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transverse  axis  were  modeled  as  a  series  connection.  The  field  gradients  and  inductances  were 
calculated  from  the  resulting  wire  patterns  using  the  Biot-Savart  law  and  Neumann  formula, 

respectively  (16). 

Fabrication 

The  fabrication  process  was  developed  and  carried  out  completely  within  our  labs  over  a  three 
month  period.  Predominantly  carried  out  by  a  single  individual,  the  estimated  total  cost  of  the 
project  was  $20,000  dollars  (CAD).  The  general  procedure  has  been  developed  during  the  con¬ 
struction  of  several  previous  gradient  coils  and  was  customized  as  required  for  this  project.  Of 
special  concern  was  the  production  of  a  very  robust  gradient  coil  capable  of  operation  at  high  cur¬ 
rents  within  a  typical  clinical  scanner.  For  this  reason  special  attention  was  paid  to  the  potting 
procedures  as  described  below. 

The  wire  used  was  composed  of  1 19  strands  of  36  AWG  copper  wire  arranged  in  a  Litz  braid¬ 
ing  to  produce  a  rKtangular  cross  section  wire  0.8  x  3.5  nun^  in  size  with  an  equivalent  gauge  of 
16  AWG  (New  England  Electric  Wire  Corporation,  New  England,  USA).  Individual  conducting 
elemente  were  coated  with  insulation  heat  rated  to  150°  C.  Litz  wire  was  used  in  place  of  solid 
copper  to  reduce  resistance  at  frequencies  of  several  kHz,  expected  to  be  the  dominant  operating 
frequency  of  the  coil  for  high  speed  imaging  sequences.  The  DC  resistance  of  the  wire  was 
BmO/m.  The  multistrand  formation  of  the  Litz  wire  also  improved  its  workability,  reducing  wire 
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memory  and  allowing  a  minimum  bending  radii  of  several  mm.  In  comparison,  winding  the  trans¬ 
verse  axis  patterns  with  equal  area  solid  copper  wire  was  found  to  be  exceedingly  difficult. 

The  wire  patterns  for  each  layer  were  CNC  milled  into  the  outer  surfaces  of  teflon  cylinders 
using  a  3 -axis  cylindrical  mill  developed  and  constructed  within  our  labs  for  the  specific  purpose 
of  gradient  coil  fabrication.  A  complete  pattern  for  a  single  transverse  layer  could  be  cut  in 
approximately  6  hours  in  a  single  pass.  Grooves  were  milled  1.75mm  deep  (half  the  height  of  the 
wire)  and  precisely  0.8mm  wide  such  that  when  the  wire  was  pushed  into  the  grooves,  it  would  be 
held  in  the  plastic  former  while  extending  1.75mm  out  of  the  cylinder  wall. 

The  winding  of  the  layers  proceeded  relatively  slowly,  typically  taking  a  single  person  an 
average  of  18  hours  to  finalize  each  layer.  After  winding  each  layer,  a  light  coating  of  epoxy  was 
applied  to  help  hold  the  wire  in  place  and  to  fix  a  type-E  thermocouple  over  the  wire  at  the  loca¬ 
tion  of  maximum  wire  density.  These  thermocouples  allow  active  temperature  monitoring  of  the 

hottest  areas  of  the  coil,  thereby  permitting  thermal  characterization  of  the  coil  as  discussed 

\ 

below,  as  well  as  the  continual  monitoring  of  the  coil  during  operation  to  prevent  accidental  dam¬ 
age.  'Pie  pattern  for  each  layer  was  wound  with  one  continuous  length  of  wire  to  avoid  the  pres¬ 
ence  of  solder  connections  within  the  wall  of  the  coil,  the  most  probable  site  for  thermal  failure  in 
our  experience.  The  five  wound  layers  were  then  potted  together  into  one  single  cylinder.  A  1.5 
cm  thick  layer  containing  7mm  diameter  copper  tubing  carrying  forced  water  cooling  at  flow  rates 
of  up  to  5l/min  was  potted  around  the  exterior  of  the  coil  to  form  the  outermost  layer.  The  epoxy 
used  for  the  potting  (Durapot  865,  Cotronics  Corporation,  New  York,  USA)  was  specially  chosen 
to  have  a  very  high  thermal  conductivity  (2.8W/m/K)  so  that  the  cylinder  wall  could  effectively 
transfer  heat  from  the  wire  layers  to  the  outer  cooling  tubes.  The  epoxy  was  heated  before  potting 
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to  reduce  its  viscosity  and  the  entire  coil  structure  was  mechanically  agitated  during  the  first  few 
hours  of  curing  to  assist  the  penetration  of  epoxy  and  prevent  air  pocket  formation.  The  potting 
was  accomplished  in  stages,  one  for  each  layer,  with  each  stage  requiring  a  minimum  of  24  hours 

to  allow  the  epoxy  to  cure. 

The  aperture  designed  for  the  side  of  the  coil  was  maintained  during  the  pottmg  procedure  and 
then  honed  as  a  final  step  to  assure  an  aperture  clear  bore  diameter  of  1cm.  Figure  4  shows  a  pho¬ 
tograph  of  the  completed  coil  before  the  addition  of  the  outermost  cooling  layer.  Connections 
between  layers  were  made  via  heavy-duty  copper  jumpers  at  the  back  of  the  coil.  Thermocouple 
lead  connectors  were  potted  into  the  back  wall  of  the  coil  for  stability.  The  inner  diameter  is  5cm, 
outer  diameter  is  15cm,  and  total  coil  length  (including  cable  connection  posts  and  cooling  con¬ 
nectors)  is  25cm.  The  completed  coil  weighs  approximately  23lbs. 

Heating 

Using  thermocouples  incorporated  into  the  coil  during  construction,  internal  temperature  mea- 

■■  X , 

surements  were  madeWer  time  while  passing  20A  DC  through  each  axis  in  turn.  Heatmg  was 
halted  when  internal  coil  temperatures  increased  by  30°C .  The  cooling  rate  of  the  coil  was  mea¬ 
sured  by  recording  the  internal  temperature  as  a  function  of  time  after  the  cessation  of  current. 
The  experiment  was  repeated  for  each  axis  for  water  flow  rates  of  01/min  and  41/mm.  Both 
and  T  were  obtained  by  applying  a  least-squares  fit  to  the  experimental  curves  for  each  axis  (equa¬ 
tion  [5]),  during  both  heating  and  cooling  phases. 
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Electrical  Measurements  and  Imaging 

Bench  top  measurements  of  both  resistance  and  inductance  were  made  at  frequencies  of  IkHz 
and  lOkHz  using  a  Hewlett-Packard  4262A  LCR  meter. 

The  completed  coil  was  directly  interfaced  with  both  a  Varian-Siemens  Unity  INOVA  4T 
whole  body  scanner  (peak  current:  250A,  peak  voltage:  600V)  and  a  1.5T  General  Electric  Signa 
Horizon  with  Echospeed  gradient  amplifiers  (peak  current:  200 A,  peak  voltage:  1200V,  General 
ElectricMedical  Systems,  Milwaukee,  WI).  No  eddy  current  compensation  was  applied  for  any 
of  the  sequences  run.  The  coil  was  well  fixated  within  the  magnet  bore  in  each  circumstance. 
Individual  solenoidal  RF  coils  (18mm  and  34mm  in  diameter)  were  constructed  and  tuned  to  res¬ 
onate  within  the  gradient  coil  bore  for  both  field  strengths. 

The  gradient  efficiency  was  measured  by  imaging  a  precision  grid  phantom  (square  in  plane 
grid  elements  of  size  2x2  mm^)  using  a  known  current.  After  reconstruction,  the  image  measure¬ 
ments  were  used  in  conjunction  with  the  known  phantom  dimensions  to  adjust  the  gradient  effi¬ 
ciency  until  the  grid  elements  in  the  centre  of  the  image  were  of  the  correct  dimension. 

Results 

"‘'■Sv-'- 

Optimized  Layer  Design  and  Interleaving 

Figure  5  shows  a  contour  plot  of  rj  in  mT/m/A  over  the  parameter  axes  of  ROU  and  coil  length 
for  a  fixed  winding  density  of  8cm'^  and  radius  of  2.775cm.  The  ROU  for  the  purposes  of  this 
plot  is  taken  to  be  the  longitudinal  length  of  the  region  within  which  the  absolute  value  of  the  gra¬ 
dient  inhomogeneity  is  less  than  30%;  however,  the  relative  features  of  the  plot  are  unchanged 
with  the  use  of  more  stringent  homogeneity  specifications  defining  the  ROU.  For  aspect  ratios 
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below  approximately  1.5  (coil  length  of  approximately  8cm  on  the  plot),  the  general  trend  is  that 
the  most  efficient  designs  correspond  to  the  regions  of  smallest  ROU  and  longest  coil,  however, 
the  trend  for  aspect  ratios  above  about  1.6  (coil  length  of  approximately  9cm  on  the  plot)  is  that 
the  designs  get  more  efficient  as  both  the  coil  length  and  ROU  are  increased.  The  dependency  of 
T]  on  coil  length  (aspect  ratio)  is  essentially  monotonic  over  the  entire  plot  with  increasing  length 

resulting  in  more  efficient  designs. 

For  mouse  imaging  with  minimal  gradient  distortion  over  the  heart  region  (<lcm  in  diameter 
(10))  and  correctable  distortion  over  a  region  large  enough  to  image  the  mouse  body  if  desired 
(~3cm  in  diameter  (10))  a  region  of  better  than  1%  homogeneity  was  required  over  1cm  and  30% 
over  3cm  diameter  spherical  volumes  (DSV).  Taking  the  30%  DSV  as  the  primary  homogeneity 
constraint,  a  line  can  be  drawn  on  the  design  space  contour  plot  (figure  5)  to  indicate  the  region  of 
coils  that  can  be  considered.  If  there  were  length  constraints  for  this  coil,  a  horizontal  Ime  at  that 
maximum  length  could  also  be  added  to  fiirther  limit  the  region  under  consideration.  In  this  appli¬ 
cation  there  is  no  sucli'constraint  and  the  entire  plot  to  the  right  of  the  DSV  line  is  allowed.  The 
chosen  design  is  at  the  point  defined  by  a  coil  length  of  1 1.1cm  and  an  ROU  size  of  3.3cm  (indi¬ 
cated  by  a  star  on  figure  5).  The  higher  efficiency  region  above  and  to  the  right  of  the  selected 
design  was  excluded  because  of  increased  inductance  values  associated  with  that  area  of  the 
space.  The  optimal  layer  design  for  the  longitudinal  axis  was  obtained  in  the  same  manner.  The 
calculated  values  of  q,  L,  and  R  are  listed  in  the  first  row  of  table  1  for  the  innermost  winding 
radius  of  2.775cm.  All  other  layers  of  both  the  X  and  Y  axes  were  based  on  the  design  of  figure 
6a,  each  scaled  to  the  appropriate  radius  while  maintaining  an  approximately  constant  winding 
density  of  8cm'’.  The  calculated  performance  values  for  each  of  these  layers  are  also  shown  m 


19 


table  1 .  The  wire  patterns  for  a  single  octant  of  the  transverse  and  longitudinal  layers  are  shown  in 
figures  6a  and  6b  respectively.  The  layer  ordering  implicit  in  table  1  is  the  optimal  one  as  deter¬ 
mined  from  the  interleaving  results  described  below. 

Table  2  lists  the  calculated  results  for  the  possible  layer  orderings.  The  most  unbalanced  case 
(XXZYY)  results  in  a  range  of  efficiency  of  74%  (X  largest,  Y  smallest)  and  a  range  of  inductance 
of  332%  (Y  largest,  X  smallest).  For  the  most  balanced  ordering  (XYZYX),  the  efficiency  and 
inductance  ranges  were  reduced  to  14%  and  25%,  respectively.  In  this  application  it  was  decided 
that  the  most  balanced  combination  was  best  and  therefore  the  XYZYX  interleave  was  used. 

We  also  investigated  the  benefit  of  adding  additional  layers  to  the  designs.  All  possibilities 
were  modeled,  and  the  results  summarized  in  table  3.  Looking  at  the  most  efficiency  balanced 
example  (XYZYXYXZ),  it  is  seen  that  the  average  efficiency  of  the  transverse  axes  increases  by 
29%  while  the  average  inductance  increases  by  242%.  This  may  be  warranted  in  situations  where 
switching  speed  is  not  important  or  where  the  larger  inductance  is  needed  to  simplify  interfacing 
the  coil  with  existing- amplifiers;  however,  in  this  application  it  was  decided  that  the  increases 
were  not  useful. 

V 

Electrical  Performance 

The  measured  values  for  resistance  and  inductance  of  the  individual  layers  as  well  as  the  com¬ 
plete  axes  are  summarized  in  table  4.  The  calculated  values  are  also  included  for  comparison.  It 
should  be  noted  that  the  relatively  slow  increase  of  resistance  with  frequency  is  expected  because 
of  the  Litz  wire  construction.  The  calculated  resistance  values  are  generally  lower  than  the  exper¬ 
imental  results  because  they  assume  zero  frequency  (DC)  and  they  do  not  take  into  account  the 
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length  of  connecting  wire  within  the  coil  or  the  effects  of  linking  the  wire  to  the  brass  connection 
posts.  These  effects  are  negligible  for  inductance,  as  indicated  by  the  accuracy  of  the  inductance 
calculations  as  compared  to  the  measurements. 

Heating 

The  results  of  temperature  measurements  as  a  function  of  time  for  the  x  axis  during  periods  of 
heating  and  subsequent  cooling  are  shown  in  figure  7.  Superimposed  on  the  data  points  are  the 
best-fit  curves  based  on  equation  [5].  The  heating  curves  for  both  the  y  and  z  axes  are  similar. 
The  fitted  values  of  and  x  for  all  measurements  are  summarized  in  table  5.  Table  5  also 
includes  the  corresponding  values  of  as  defined  in  equation  [6].  The  data  indicates  that  the 
steady  state  heating  parameter,  k^,  is  decreased  by  an  average  of  46%  with  the  addition  of  water 
cooling  (46%,  21%  and  70%  for  the  X,Y,  and  Z  axes,  respectively).  The  heating  time  constant 
(Xh)  is  always  found  to  decrease  with  cooling,  whereas  the  cooling  time  constant  (Xc)  always 
increases  (although  only  very  slightly  for  the  z  axis).  The  most  important  effect  of  coolmg  is  evi¬ 
dent  by  noting  the  difference  in  steady  state  temperature  decreases.  These  results  suggest  that 
without  cooling  the  steady  state  temperatures  of  the  coil,  after  the  cessation  of  power,  approach  on 
average' approximately  15°C  above  ambient  temperature.  Of  course,  eventually  the  coil  would 
cool  to  room  temperature;  however,  this  full  return  to  baseline  temperatures  required  several 
hours.  In  practice,  this  would  represent  a  serious  limitation  to  coil  performance.  The  addition  of 
cooling  removed  this  limitation  and  resulted  in  cooling  times  very  close  to  heating  times  (noted  by 
the  similarity  between  x^  and  x,.  with  the  inclusion  of  cooling). 

Figure  8  shows  the  number  of  single  shot  EPI  images  that  can  be  collected  at  maximum  rate 
(as  determined  by  the  pulse  sequence)  during  each  burst  period  as  a  function  of  FOV  for  each 
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axis.  The  x  axis  allows  the  collection  of  the  largest  number  of  images  before  each  dead  time  (14% 
more  than  the  z  axis,  57%  more  than  the  y  axis);  however,  this  is  balanced  by  the  fact  that  the 
entire  burst  plus  dead  time  sequence  can  be  repeated  every  430  seconds  when  using  the  z  axis,  as 
compared  to  631  and  666  seconds  for  the  Y  and  X  axes,  respectively.  In  practice,  burst  mode  can¬ 
not  be  accomplished  repeatedly  without  water  cooling  because  of  the  extremely  long  times 
required  for  the  coil  to  cool. 

Figure  9  shows  a  plot  of  the  minimum  single  shot  EPI  repetition  time  as  a  function  of  the  read¬ 
out  FOV  for  the  three  individual  axes.  The  z  axis  allows  an  imaging  rate  approximately  two  times 
that  of  y  and  50%  more  than  x  over  the  full  range  of  FOVs,  while  the  x  axis  demonstrates  an 
approximate  30%  speed  advantage  compared  to  the  y.  A  similar  analysis  can  be  used  to  compare 
the  performance  of  each  axis  with  and  without  water  cooling  and  it  was  foimd  that  the  z  axis 
showed  a  decrease  in  minimum  TR  by  70%,  the  x  axis  by  45%,  and  the  y  axis  by  20%  with  the 
use  of  water  cooling. 

-  V 

Preliminary  Imaging  • 

Figure  10  shows  the  first  image  collected,  on  the  4T  scanner,  of  the  grid  phantom  as  described 

V 

in  the  Methods  section.  It  was  collected  with  a  gradient-recalled  echo  sequence  with  the  follow¬ 
ing  parameters:  TR  =  50ms,  TE  =  4ms,  flip  angle  =11°,  matrix  size  =  256  x  256,  FOV  =  4x4 
cm^,  slice  thickness  =  2mm,  and  number  of  averages  =16.  The  image  in  the  figure  is  cropped  to 
highlight  the  grid  phantom  of  size  26mm  square.  A  box  has  been  overlaid  on  the  image  to  indi¬ 
cate  a  region  of  1cm  square  that  corresponds  to  the  size  of  the  ROU  over  which  less  than  1%  devi¬ 
ation  was  expected.  It  can  be  seen  that  no  obvious  image  distortions  occur  over  that  area.  The 
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measured  efficiencies  at  the  coil  center  for  the  X,  Y,  and  Z  axes  were  found  to  be  1 1.8  mX/m/A, 

1 1 .0  mX/m/A,  and  1 1 .2  mX/m/A  respectively. 

Figure  1 1  shows  a  proton-density-weighted  fast  spin-echo  image  of  an  excised  rat  brain  col¬ 
lected  on  a  GE  1 .5X  Signa  scanner.  Xhe  sequence  parameters  were:  XR  =  3500ms,  effective  XE  - 
18ms,  matrix  size  =  256  x  192,  FOV  =  1.7  x  1.7  cm^,  slice  thickness  =  0.625mm,  echo  train  length 
=  4,  and  number  of  averages  =  24.  Xotal  imaging  time  was  approximately  one  hour  and  fifteen 
minutes.  Xhe  in-plane  resolution  was  65  x  87  pm^.  Xhe  smooth,  round  outer  surface  of  the  brain 
is  due  to  the  14mm  diameter  tube  into  which  the  brain  was  placed  during  imaging. 

Discussion 

Although  there  are  numerous  publications  describing  gradient  coil  designs  (17),  descnptions 
of  coils  customized  for  small  animal  imaging  are  lacking.  A  body  of  work  does  exist  for  the 
development  of  very  small,  extremely  high  strength  designs  for  use  in  NMR  diffusion  studies  m 
very  small  (typically  less  than  Icm^)  samples  (18-28).  Xhe  very  small  coils  described  for  these 
studies  are  typically  one  or  two  axis  designs  and  hence  are  not  suitable  for  microimaging  studies 
(22-24).  Quadrupolar  designs  (for  at  most  two  axes)  are  the  most  popular  in  these  applications 
(19, 20, 23-25, 29, 33, 34),  assuming  that  there  is  adequate  room  within  the  scanner  bore  for  trans¬ 
verse  coil  alignment.  Coils  in  the  transverse  orientation  (29,  33-35)  and  on-axis  (30-32,36)  have 
been  reported  that  are  large  enough  be  used  for  small  animal  studies  while  multilayer  designs  for 
single  axes  have  been  reported  (26-28)  for  diffusion  NMR  studies,  but  there  are  no  complete  three 
axis  designs  considered  that  apply  multilayer  design  methods,  nor  have  multilayer  methods  been 
applied  to  the  design  of  coils  large  enough  for  small  animal  imaging.  A  survey  of  the  literature 
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indicates  that  the  strongest  on-axis  gradient  coils  in  use  for  mouse  imaging  produce  gradients  of 
up  to  lOOOmT/m  (10).  The  trade-offs  between  coil  aspect  ratio,  imaging  ROU,  and  coil  perfor¬ 
mance  have  not  been  demonstrated  previously,  nor  have  the  special  considerations  that  must  be 
taken  into  account  to  properly  govern  coil  heating. 

This  paper  demonstrates  that  the  primary  challenges  to  attaining  the  extreme  gradient  strength 
and  speed  requirements  that  accompany  mouse  imaging  pertain  to  the  thermal  response  of  the 
coil.  To  alleviate  the  extreme  heating  problems  predicted  by  the  thermal  scaling  laws  described  in 
the  theory  section,  efforts  were  made  on  several  fronts  to  manage  the  temperature  of  the  coil  while 
maintaining  the  required  strength  and  inductance.  A  multilayer  approach  was  taken  that  allowed 
high  efficiency  designs  with  reduced  winding  density.  The  lower  winding  densities  are  vital  in 
that  they  allow  larger  cross-section  wire  to  be  used,  thereby  reducing  resistance  significantly.  We 
developed  fabrication  methods  that  allow  the  use  of  rectangular  Litz  wire.  The  rectangular  wire 
was  wound  on  end,  increasing  the  available  copper  area  by  a  factor  of  four  (for  the  particular  wire 
used  in  this  coil)  as  compared  to  square  wire  wound  to  the  same  density.  The  Litz  wire  was  cho¬ 
sen  for  its  resistance  characteristics  over  the  1  to  10  kHz  range,  expected  to  be  the  dominant  oper¬ 
ating  frequency  for  the  high  speed  image  acquisitions  intended  for  this  coil.  These  measures  all 
contributed  to  reducing  the  resistance  of  each  axis.  Going  further,  we  sought  to  reduce  the  tem¬ 
perature  increases  for  any  input  power  through  the  addition  of  water  cooling.  The  epoxy  used  for 
the  fabrication  of  the  coil  was  specially  chosen  to  have  the  highest  thermal  conductivity  available. 
This  maximized  the  transfer  of  heat  from  the  wires  to  the  water  cooling  layer.  The  potting  proce¬ 
dure  was  specially  developed  to  avoid  the  inclusion  of  any  heat  insulating  materials  that  would 
hinder  the  removal  of  heat  from  the  coil.  Thermocouples  were  built  into  the  coil  to  allow  both  the 
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characterization  of  the  heating  response  of  the  coil  under  various  circumstances,  as  well  as  to 
facilitate  continual  thermal  monitoring  during  operation.  Finally,  efforts  were  made  to  model  the 
heating  effects  of  two  very  different  modes  of  operation.  These  efforts  all  contributed  to  attaining 

Optimal  performance  with  the  final  coiL 

It  is  instructive  to  compare  the  thermal  response  as  predicted  by  the  imtial  scaling  laws  to  the 
measured  response  of  the  completed  coil.  From  figure  1,  we  saw  that  if  the  hypothetical  body  coil 
were  simply  scaled  down  to  mouse  size,  a  winding  density  of  16cm*‘  would  have  been  required  to 
achieve  the  target  efiiciency.  The  corresponding  steady-state  temperature  increase  at  100 A  was 
seen  to  be  approximately  twenty-five  thousand  degrees,  clearly  not  feasible  and  underlming  the 
need  for  a  different  approach.  From  the  heating  measurements  made  and  summarized  m  table  5, 
we  saw  that  for  an  input  rms  current  of  20A,  the  average  steady-state  temperature  increase  of  the 
axes  was  32«C .  Scaling  this  result  to  an  input  rms  current  of  lOOA,  the  temperature  increase  of 
the  axes  is  predicted  to  be  approximately  800°C.  Although  still  a  much  larger  temperature 
increase  than  tolerable^  it  is  more  than  thirty  times  lower  than  that  predicted  by  the  umform  scal¬ 
ing  law.  This  substantial  improvement  is  due  to  the  many  steps  taken  to  limit  the  heatmg  of  this 

mouse  coil. 

The  results  indicate  that  this  coil  design  can  be  successfully  operated  in  both  burst  and  contin¬ 
uous  modes.  In  different  applications,  both  will  be  necessary  and  it  is  important  to  be  able  to  opti¬ 
mize  the  coil  duty  cycle  for  either  case.  Burst  mode  will  allow  the  collection  of  hundreds  of 
successive  <lcm  FOV  single  shot  64  x  64  EPI  images  at  rates  of  approximately  60  images  per 
second.  Continuous  mode  will  allow  steady  collection  of  <lcm  FOV  single  shot  64  x  64  EPI 
images  at  rates  of  approximately  one  per  second.  These  are  of  course  the  extreme  cases  with 
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respect  to  FOV.  Practical  operation  will  undoubtably  involve  some  compromise  between  FOV 
and  imaging  rates;  the  theory  and  results  described  in  this  paper  can  be  applied  to  optimize  the 
coil  operation  in  any  circumstance. 

There  are  several  ways  that  the  coil  design  presented  here  can  be  improved.  Firstly,  minimum 
power  instead  of  minimum  inductance  methods  could  have  been  applied  and  this  may  have 
reduced  the  axis  resistance  at  the  cost  of  increased  inductance  (17).  The  reason  that  this  was  not 
done  was  that  we  have  not  implemented  minimum  power  mediods  in  conjunction  with  current 
constraints,  as  we  have  with  CCMI.  The  use  of  current  constraints  is  important  to  control  the  final 
size  of  the  coil  and  was  necessary  to  achieve  the  side-hole  aperture  specified.  Secondly,  more  ele¬ 
gant  minimum  power  3D  methods  for  multilayer  gradient  coils  have  been  reported  recently 
(26,28).  It  would  be  possible  to  utilize  those  techniques  to  obtain  an  equivalent  strength  coil  with 
lower  resistance  in  place  of  simply  scaling  the  size  of  a  single  design  to  the  different  layers. 
Again,  this  is  expected  to  be  at  the  cost  of  increased  inductance.  We  focused  primarily  on 
improved  fabrication  i^ethods  to  minimize  the  thermal  problems  in  this  coil;  however,  because 
temperature  response  remains  the  dominant  limiting  factor  in  this  coil,  additional  design  optimi- 

I 

zation  efforts  are  warranted. 

Conclusion 

This  paper  has  demonstrated  the  most  important  issues  in  the  design  and  operation  of  a  very 
high  strength  gradient  coil  for  small  animal  microscopy.  With  current  generation  amplifier  hard¬ 
ware,  this  coil  allows  the  collection  of  single  shot  EPI  images  at  6mm  FOV  and  94pm  isotropic 
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voxels  at  imaging  rates  exceeding  SOs'^.  This  coil  should  allow  improvements  in  in  vivo  MR 
microscopy  for  a  diverse  range  of  applications. 
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Figure  Captions: 

FIG.  1.  Plots  of  (a)  gradient  coil  efficiency,  (b)  inductance,  and  (c)  average  temperature  rise,  as 
functions  of  radius  for  five  different  fixed  values  of  maximum  winding  density.  The  hypothetical 
body  coil  design  upon  which  all  the  scalings  are  based  is  indicated  by  the  cross  and  the  effective 
radius  of  the  proposed  mouse  coil  is  represented  by  the  vertical  line.  The  position  of  the  target 
efficiency  (a)  and  maximum  allowed  inductance  (b)  are  indicated  by  the  horizontal  lines  on  the 
respective  plots. 

FIG.  2.  Schematic  diagrams  of  the  temperature  variation  during  (a)  burst  mode  and  (b)  continu¬ 
ous  mode  operation.  Solid  vertical  bars  at  the  bottom  of  each  diagram  represent  the  acquisition  of 
a  single  EPI  image. 

FIG.  3.  Schematic  diagrams  of  five  layer,  three  axis  coils  interleaved  in  (a)  unbalanced  and  (b) 
balanced  configurations. 

FIG.  4.  Photograph  of  gradient  coil  prior  to  addition  of  outer  cooling  layer.  Electrical  jumpers 
connect  layers  within  the  X  and  Y  axes  in  series.  Six  free  brass  electrical  cormections  (two  per 
axis)  are  visible  at  the  top  of  the  coil.  Four  thermocouple  connectors  are  visible  at  the  front  of  the 
top  surface.  The  1cm  aperture  through  the  side  is  visible  halfway  up  the  front  of  the  coil. 

FIG.  5.  Contour  plot  of  the  efficiency  [mT/m/A]  of  the  innermost  transverse  layer  (radius  = 
2.775cm)  over  axes  of  coil  length  and  30%  ROU  size.  The  position  of  the  chosen  transverse  layer 
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design  is  indicated  by  the  star.  The  vertical  line  marks  the  minimal  allowed  ROU  size  for  mouse 
imaging. 

FIG.  6.  Single  octant  wire  patterns  for  the  (a)  transverse  and  (b)  longitudinal  axes.  Arrowheads 
indicate  relative  direction  of  current  within  each  pattern.  A  region  of  current  reversal  is  present  on 
the  longitudinal  (z)  axis  only.  The  pattern  in  (a)  is  used  for  all  transverse  layers  within  the  coil, 
shifted  by  ninety  degrees  for  the  y  axis  layers.  The  maximum  winding  density  on  the  transverse 
and  longitudinal  patterns  are  8cm  ^  and  10cm  ,  respectively. 

FIG.  7.  Measured  temperatures  at  the  region  of  highest  wire  density  on  the  X-axis  over  time 
while  passing  20A  DC  through  the  X-axis  alone.  Heating  was  stopped  when  the  maximum  tem¬ 
perature  increased  by  SO^C .  The  solid  triangles  represent  the  response  with  no  cooling  applied, 
wheras  the  solid  diamonds  represent  the  response  with  41/min  water  cooling.  The  solid  lines  are 
the  results  of  fitting  equation  (5)  to  the  measured  data.  The  response  of  the  Y  and  Z  axes  were 

similar. 

FIG.  8.  The  number  of  single  shot  EPI  images  (Np  =  64,  =  64,  SW  =  500kHz,  V  =  400V,  ATj, 

=  ATni  =  10°C ,  water  cooling  =  41/min,  and  tip  Lp  R,  and  kj,  equal  to  the  values  for  each  axis  as 
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given  in  tables  2  and  4)  that  can  be  collected  by  each  axis  during  a  single  burst  period,  as  a  func¬ 
tion  of  field  of  view. 

FIG.  9.  The  minimum  TR  for  a  single  shot  EPI  image  (Np  =  64,  =  64,  SW  =  500kHz,  V  = 

400V,  ATc  =  10°C ,  water  cooling  =  41/min,  and  r|p  Lp  R,  and  kj^  equal  to  the  values  for  each  axis 
as  given  in  tables  2  and  4)  for  continuous  mode  scaiming,  as  a  function  of  field  of  view. 

FIG.  10.  Gradient  echo  image  of  a  grid  phantom  collected  with  the  mouse  gradient  coil  at  4T. 
Imaging  parameters  were:  TR  /  TE  /  flip  angle  =  50ms  /  4ms  / 1 1  ° ,  number  of  averages  =  16,  FOV 
=  4  X  4cm^,  matrix  =  256  x  256,  resolution  =  156  x  156  x  2000pm^,  collection  time  =  3min  24sec. 
The  image  has  been  cropped  to  highlight  the  phantom  of  size  26  x  26mm^.  Dark  squares  within 
the  grid  are  2  x  2mm^  and  are  separated  by  2mm.  The  dashed  box  represents  an  area  10  x  lOmm^ 
within  which  the  gradient  is  uniform  to  better  than  1%. 

FIG.  1 1 .  Fast  spin  echo  image  of  an  excised  rat  brain  collected  at  1 .5T.  Imaging  parameters  were: 
TR  /  TE  =  3500ms  /  18ms,  echo  train  length  =  4,  number  of  averages  =  24,  FOV  =  17  x  17mm^, 
matrix  =  256  x  192,  resolution  =  65  x  87  x  625iJ,m^,  collection  time  =  Ihr  15min.  '  The  smooth, 
round  duter  surface  of  the  brain  is  due  to  the  14mm  diameter  tube  into  which  the  brain  was  placed 
during  imaging. 
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radius  [cm] 


aw[cm'^] 


8.3 


r\  [mT/m/A] 


.8 


R  [mQ] 

12 

67 

23 

104 

97 

200 

60 

201 

98 

275 

Table  2;  Calculated  Results  for  Interleaving  Combinations 
“  n  [mT/m/A] 


XXYYZ 


XYXYZ 


XYYXZ 


XXYZY 


XYXZY 


XYYZX 


XZXYY 


XZYXY 


XZYYX 


ZXXYY 


Typnc/ 

ZXYYX 


60 

161 

231 

64 

114 

231 

84 

102 

231 

60 

182 

154 

64 

141 

154 

121 

102 

154 

84 

141 

97 

64 

259 

56 

84 

182 

56 

121 

161 

56  1 

102 

259 

28 

114 

182 

28 

141 

161 

28 

31 


Table  3:  Calculated  Values  for  Increased  Number  of  Layers 


GX 

[mT/m/A] 

GY 

GZ 

GX 

L[nH] 

GY 

GZ 

XYZYX 

11.7 

10.3 

11.0 

121 

114 

97 

XYZYXXYZ 

15.1 

13.4 

18.9 

419 

396 

1064 

XYZYXXZY 

15.1 

13.1 

19.7 

419 

444 

891 

XYZYXZXY 

14.7 

13.1 

20.1 

445 

444 

701 

XYZYXYZX 

14.4 

13.7 

19.7 

488 

363 

891 

XYZYXZYX 

14.4 

13.4 

20.1 

488 

396 

701 

Table  4:  Calculated  and  Measured  Electrical  Values 


calc. 

1kHz 

10kHz 

calc. 

R  [mQ] 
1kHz 

10kHz 

layer  1 

12 

11.9 

11.8 

67 

90 

96 

layer  2 

23 

23.0 

23.0 

104 

124 

131 

layer  3 

97 

98.0 

97.5 

200 

238 

260 

layer  4 

60 

61.0 

61.0 

201 

262 

273 

layer  5 

98 

101.0 

101.0 

275 

323 

345 

Xaxis 

121 

126.3 

126.1 

342 

412 

429 

Y  axis 

114 

118.4 

118.2 

305 

380 

400 

Zaxis 

97 

100.0 

99.7 

200 

238 

254 

Table  5:  Fit  Results  for  Heating  Measurements 


Axis  Flow  Phase 

T[S] 

kh[”CA-2] 

heating 

49 

1053 

0.297 

ui/min 

w  cooling 

-11 

236 

heating 

34 

666 

0.204 

41/min 

cooling 

-31 

665 

heating 

51 

786 

0.332 

ui/min 

Y  cooling 

-17 

317 

heating 

40 

631 

0.274 

41/min 

cooling 

-31 

649 

heating 

37 

768 

0.388 

ui/min  - 

-j  cooling 

-17 

253 

heating 

21 

430 

0.228 

41/min 

cooling 

-22 

256 

[V/UJ/iiu]  U 


[HUill 


0.001 
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Appendix  II 


The  IVIM  method  and  the  analysis  program  developed  for  this  grant  were  tested 
with  two  substances  -  water  and  acetone  -  whose  diffusion  coefficients  are  well 
documented.  The  theoretical  and  experimentally  determined  results  are  shown  below  for 
diffusion  encoding  along  six  different  directions: 

Water: 


The  theoretical  value  of  the  diffusion  coefficient  of  water  is  between  2.12  x  10'^ 
and  2.63  x  10'^  mm^/s  at  21.7  C°. 


Diffusion  Encoding  Direction 

Diffusion  Coefficient  tmm^/s) 

xy 

2.03  +  0.02  X  10'^ 

-xy 

2.10  ±  0.09  X  10'" 

xz 

2.08  ±  0.10  X  10'^ 

-xz 

1.97  ±  0.06  X  10'^ 

yz 

2.07  ±  0.04  X  10'^ 

-yz 

2.07  ±0.1  lx  10'^ 

Acetone: 


•  “3 

The  theoretical  value  of  the  diffusion  coefficient  of  acetone  is  between  3.88  x  10’ 
and  4.43  x  10'^  mm^/s  at  21.7  C°. 


Diffusion  Encoding  Direction 

Diffusion  Coefficient  tmmVs) 

xy 

4.47  ±  0.04  X  10'^ 

-xy 

4.59  ±  0.05  X  10'" 

XZ 

4.62  ±  0.03  X  10'^ 

-xz 

4.30  ±  0.05  X  10'^ 

yz 

4.47  ±  0.06  X  10'^ 

-yz 

4.55  ±  0.06  X  10’^ 

Reference:  Bammer  R,  et  al.  Diffusion-weighted  imaging  with  navigated  interleaved 
echo-planar  imaging  and  a  conventional  gradient  system.  Radiology  211:799  (1999). 


Diffusion  Coefficients  (xlO  mm  /s) 


Shown  below  are  the  bi-exponential  diffusion  coefficients  before  (Da,  Db)  and  after  (Dx, 
Dy)  cessation  of  blood  flow  (euthanasia)  in  non-necrotic  regions  of  interest  in  low  and 
high  dose  tumours. 


Note  that  the  fast  decay  component  (large  diffusion  coefficient),  Db,  disappears  post 
cessation  of  blood  flow.  The  presence  of  a  bi-exponential  decay  post  cessation  of  blood 
flow  may  be  attributed  to  different  rates  of  diffusion  in  the  intra-  and  extracellular  spaces. 
We  are  currently  endeavoring  to  assign  these  two  diffusion  components. 
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Practical  Design  of  a  High-Strength  Breast  Gradient  Coil 

Cynthia  F.  Maier,  Hristo  N.  Nikolov,  Kenneth  C.  Chu,  Blaine  A.  Chronik,  Brian  K.  Rutt 


A  high-strength  three-axis  local  gradient  coil  set  was  con¬ 
structed  for  MRI  of  the  breast.  Gradient  fields  with  good  uni¬ 
formity  (<10%  deviation  from  the  desired  gradient)  over  most 
of  the  volume  required  for  breast  imaging  were  generated 
with  efficiencies  of  up  to  3.3  mT/m/A.  The  coils  will  allow 
diffusion  breast  imaging  in  clinically  acceptable  examination 
times.  The  electrical  design,  water  cooling  system,  and  fabri¬ 
cation  techniques  are  described.  Preliminary  tests  of  the  coil 
included  images  of  a  grid  phantom  and  diffusion  measure¬ 
ments  in  a  short-Tg  agarose  gel  phantom. 

Key  words:  breast  imaging;  diffusion  imaging;  gradient  coils. 

INTRODUCTION 

In  this  paper,  we  present  a  specialized  gradient  coil  design 
for  breast  imaging  and  demonstrate  the  utility  of  this  coil  for 
performing  MR  diffusion  imaging  with  short-duration, 
high-strength  gradient  pulses.  Sensitization  of  an  MR  pulse 
sequence  to  motion  can  be  achieved  using  conventionally 
available  gradient  hardware;  however,  measurement  of  wa¬ 
ter  diffusion  (or  slow  flow  in  microvessels)  with  this  hard¬ 
ware  requires  long  echo  times  to  achieve  large  gradient 
moments  and  therefore  suffers  from  significant  signal  loss 
due  to  T2  relaxation  and  system  nonidealities.  In  addition, 
increasing  the  time  between  motion-sensitizing  gradient 
pulses  allows  more  macroscopic  motion  to  occur  in  this 
interval,  resulting  in  pixel  misregistration,  motion-induced 
ghost  artifacts,  and  diffusion/flow  measurement  artifacts. 
Higher-strength  gradients  will  allow  adequate  sensitization 
to  small  motions  (i.e.,  large  gradient  moments)  with  very 
short,  high-amplitude  gradient  pulses  and  minimal  time 
delay  between  these  pulses. 

METHODS 

The  gradient  coil  assembly  was  designed  in  a  cylindrical 
geometry,  with  each  gradient  coil  built  onto  a  separate 
cylindrical  former.  The  three-axis  coil  set  was  composed 
of  concentrically  nested  coils,  oriented  vertically  such 
that  their  common  axis  was  transverse  to  the  main  field 
as  shown  in  Fig.  1.  Space  restrictions  within  the  magnet 
bore  require  that  the  coils  be  as  short  as  possible  in  the 
vertical  direction  to  allow  a  patient  to  lie  comfortably 
above  the  top  edge  of  the  cylinders.  Unfortunately,  de¬ 
creasing  the  height  of  the  coils  (relative  to  their  diame- 
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ters)  makes  the  creation  of  uniform  gradients  within  the 
coils  much  more  difficult.  To  satisfy  the  competing  de¬ 
mands  of  good  gradient  uniformity  and  height  restric¬ 
tions,  we  used  an  aspect  ratio  of  approximately  1:1  for 
our  cylindrical  coils.  The  height  of  the  complete  coil 
assembly  was  16  cm,  leaving  34-cm  clearance  in  the 
vertical  direction  in  our  General  Electric  (GE)  Signa  Ho¬ 
rizon  scanner  (General  Electric  Medical  Systems,  Mil¬ 
waukee,  WI),  and  the  largest  coil  diameter  was  21.6  cm. 

To  generate  the  winding  densities  required  for  high  gra¬ 
dient  strengths,  relatively  small-diameter  wire  (AWG  no. 
20)  was  used  to  construct  the  coils.  The  resistivity  of  this 
wire  is  high,  requiring  an  active  cooling  system  to  remove 
heat  and  maintain  temperatures  inside  the  coils  within  a 
range  comfortable  for  the  patient.  Water-cooling  jackets 
with  very  thin  walls  were  sandwiched  between  the  gradient 
coils,  so  that  resistive  heat  could  be  removed  efficiently 
while  maintaining  the  compact  construction  of  the  coil 
assembly.  A  dedicated  RF  coil  and  shield  were  mounted 
inside  the  innermost  gradient  coil  with  an  approximately 
1.3-cm  space  between  the  RF  coil  and  closest  gradient  coil. 
The  entire  assembly  is  depicted  in  cross-section  in  Fig.  2. 

Electrical  Design 

The  highly  efficient  quadrupolar  current  distribution  (1-3) 
was  used  to  provide  gradients  in  along  both  the  x  and  z 
directions  (G^and  GJ.  This  current  density  varies  azimuth- 
ally  as  cos (20)  and  is  entirely  directed  along  the  axial  di¬ 
rection  on  the  surface  of  the  cylindrical  former.  Discrete 
wire  positions  approximating  the  continuous  quadrupolar 
current  density  were  calculated  by  sampling  the  current 
density  cos(20)  at  0^-  such  that  cos(20j-+i^)  —  cos(20j  was 
constant.  A  schematic  diagram  showing  this  ideal  wire 
distribution  is  given  in  Fig.  3a.  The  current  pattern  in  Fig. 
3a  is  nonphysical  in  that  the  current  is  discontinuous  at  the 
ends  of  the  coil,  and  therefore,  connecting  paths  between 
the  wires  (“return  wires”)  must  be  provided.  In  most  imag¬ 
ing  applications  in  which  local  gradient  coils  have  been 
used  previously,  the  volume  to  be  imaged  has  been  located 
at  the  center  of  the  coil,  and  placing  return  paths  near  the 
ends  of  the  cylinder  had  a  minimal  effect  on  the  desired 
field  within  the  object.  For  breast  imaging  applications,  the 
required  location  for  the  region  of  good  gradient  uniformity 
is  near  the  top  edge  of  the  coils  and  a  conventional  return 
wire  arrangement  would  corrupt  the  gradient  field  unifor¬ 
mity  within  this  volume.  We  chose  to  place  the  return  wire 
components  instead  on  horizontal  plates  at  the  tops  and 
bottoms  of  the  coils.  Close  to  the  cylinders  on  the  plates,  the 
return  paths  were  extended  along  the  z  direction  for  at  least 
5  cm  (no  is  generated  by  current  in  the  z  direction).  The 
bottom  set  of  return  wires  for  the  G^  coil  were  stnmg  across 
the  bottom  of  the  cylinder  in  the  z  direction.  The  quadru- 
pole  current  distribution  for  G^  is  rotated  45°  with  respect 
to  G^,  and  the  retmn  wires  for  G^  could  not  be  strung  in  the 
z  direction;  therefore,  these  wires  were  located  on  a  hori¬ 
zontal  plate. 
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FIG.  1.  Schematic  showing  orientation  of  iocal  breast  gradient 
coil  in  magnet. 

The  design  of  the  Gy  gradient  makes  use  of  the  fact  that 
V  X  5  =  0  throughout  any  region  in  which  the  current 
density  is  0.  This  means  that  a  transverse  gradient  coil 
with  its  axis  along  z  (the  conventional  orientation)  can  be 
turned  90°  about  the  x  axis  so  that  its  axis  is  along  the  y 
direction  and  used  to  produce  a  gradient  in  along  the 
y  direction.  We  designed  a  transverse  gradient  coil  in  the 
conventional  orientation  using  the  target  field  method  (4, 
5)  and  then  rotated  the  coil  to  produce  Gy  Although  a 
wide  range  of  breast  sizes  exist,  we  defined  an  average 
volume  of  interest  (VOI)  for  breast  imaging  (by  compari¬ 
son  with  commercially  available  breast  RF  coils)  as  a 
cylindrical  volume  with  length  equal  to  60%  of  the  coil 
height  (11.5  cm),  a  diameter  equal  to  90%  of  the  RF  coil 
diameter  (13.5  cm),  and  location  flush  with  the  top  sur¬ 
face  of  the  gradient  coil  assembly.  The  target  field  coil 
was  designed  to  have  its  region  of  good  gradient  unifor¬ 
mity  centered  on  this  VOI.  This  design  consisted  of  only 
two  sets  of  closed  current  loops  (Fig.  3b).  For  the  target- 
field  coil  in  its  conventional  orientation  with  its  axis 
along  2,  the  desired  gradient  is  created  almost  entirely  by 
the  current  in  half  of  the  wire  loops  only  (the  “primary*’ 
component  of  the  loop  as  indicated  in  Fig.  3b),  and  return 
wires  are  an  intrinsic  part  of  the  design  (labeled  “return” 
in  Fig.  3b).  For  this  coil,  only  one  set  of  “return  paths”  is 
required  for  current  continuity,  and  in  the  vertical  orien- 
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10.82  cm 

FIG.  2.  Cross-section  of  complete  gradient  coll  assembly  show¬ 
ing  location  of  RF  coil  and  water-cooling  modules. 


Primary 

Wires 


Return 

Wires 


FIG.  3.  (a)  Quadrupolar  current  density  with  56  windings,  (b)  Tar- 
get-field-deslgn  transverse  coil.  The  “primary  wire"  component  is 
almost  entirely  responsible  for  generation  of  the  transverse  gradi¬ 
ent  for  the  conventional  coil  orientation  with  its  axis  along  z. 


tation,  these  wires  are  located  at  the  bottom  of  the  cylin¬ 
der.  In  our  coil  design,  such  a  large  number  of  loops  was 
required  to  generate  the  desired  gradient  efficiency  that 
these  return  wires  would  have  occupied  a  large  fraction 
of  the  space  available  on  the  surface  of  the  cylinder.  Due 
to  the  strict  limitations  on  coil  height,  we  believed  that 
this  space  could  be  used  more  efficiently  if  only  the 
primary  component  of  the  coil  design  was  present  on  the 
cylinder  itself  and  the  return  paths  were  located  on  a 
horizontal  plate  as  for  the  quadrupole  coils.  A  complica¬ 
tion  that  arises  from  truncating  the  wire  pattern  in  the 
vertical  dimension  in  this  way  is  that  the  reciprocity 
relations  from  Maxwell’s  equations  no  longer  hold  for  a 
nonphysical  current  distribution  (i.e.,  V  X  B  =5^  0  for  a 
current  distribution  with  nonzero  divergence).  This 
means  that  the  gradient  field  calculated  for  the  truncated 
design  is  not  identical  in  the  conventional  and  transverse 
orientations  unless  the  return  wires  are  included  in  the 
calculation.  We  shall  see  in  the  Results  section  that  the 
truncated  coil  in  the  transverse  orientation  has  both  an 
increased  gradient  efficiency  compared  to  the  full  design 
for  the  same  coil  height  and,  additionally,  the  region  of 
acceptable  gradient  uniformity  (<20%  deviation  from 
uniformity)  is  significantly  larger  for  the  truncated  coil. 

One  undesirable  side  effect  that  arises  from  locating 
the  return  wires  on  horizontal  surfaces  is  that  the  torque 
on  the  coils  from  interaction  with  the  main  static  field  is 
greatly  increased.  This  effect  is  significant.  For  our  1.5-T 
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magnet,  current  is  only  supplied  when  the  coils  are  po¬ 
sitioned  near  the  central,  highly  homogeneous  region  of 
the  magnet  (typical  clinical  imaging  systems  have  a  50- 
cm- diameter  spherical  volume  of  better  than  10  ppm  Bq 
uniformity)  and  the  upward  (downward)  force  on  the  +z 
(-z)  half  of  the  horizontal  plate  for  the  coil  is  approx¬ 
imately  15  N/A  (3,5  Ib/A)  increasing  linearly  with  ap¬ 
plied  current.  As  the  gradients  are  switched,  the  effect  of 
this  force  is  to  create  an  alternating  torque  of  approxi¬ 
mately  3  Nm/A  in  size.  For  many  pulse  sequences,  the 
direction  of  the  current  is  switched  rapidly,  resulting  in 
a  torque  that  alternates  direction  rapidly  and  causes  sig¬ 
nificant  vibrational  stress  on  the  coils.  For  our  prelimi¬ 
nary  experiments,  a  temporary  coil  support  was  con¬ 
structed  from  wood-  and  glass-reinforced  phenolic.  The 
horizontal  plates  were  rigidly  attached  to  the  support, 
and  the  support  structure  itself  served  to  mount  the  coil 
assembly  inside  the  bore  of  the  imager,  thereby  prevent¬ 
ing  any  bulk  coil  motions  that  could  arise  from  the  mag¬ 
netic  forces  on  the  wires.  The  wires  themselves  were 
potted  securely  in  the  milled  grooves  using  epoxy  and 
were  prevented  from  vibrating  loose  from  these  grooves 
by  the  next  adjacent  layer  of  cooling  channels,  since  all 
layers  were  tightly  sandwiched  together. 

Col!  Construction 

To  produce  high  gradient  efficiencies,  it  is  advantageous 
to  use  the  maximum  number  of  wires  possible.  This 
number  is  constrained  in  an  obvious  way  by  the  wire 
diameter,  minimum  spacing  between  wires,  and  circum¬ 
ference  of  the  cylindrical  former.  To  obtain  the  desired 
gradient  efficiencies  for  these  coils,  it  was  necessary  to 
use  relatively  small  wires  (AWG  no.  20)  with  a  diameter 
of  ~1  mm  and  an  interwire  spacing  of  0.5  mm.  The 
(outermost)  coil  was  wound  with  a  single  length  of  wire 
such  that  216  axially  directed  segments  generated  the 
required  current  density,  and  for  the  (innermost)  coil, 
200  axially  directed  segments  were  wound.  For  the  mid¬ 
dle  coil  (G^,  constraining  the  closest  separation  of  the 
wires  to  be  0.5  mm  allowed  48  loops  per  side.  For  all 
three  coils,  the  wire  paths  were  milled  (numerically  con¬ 
trolled)  into  the  surface  of  polycarbonate/acrylic  formers. 
Wire  was  then  wound  into  these  grooves  and  epoxied  in 
place.  Figure  4a  shows  the  constructed  G^  coil,  complete 
with  horizontal  plates  after  winding  with  (dielectric)  in¬ 
sulated  copper  wire.  The  upper  plates  were  curved 
downward  at  the  superior  and  inferior  ends  to  provide  a 
more  comfortable  structure  for  patient  support.  Figure  4b 
shows  the  middle  coil  (G^  before  assembly  of  the  com¬ 
plete  gradient  coil  unit. 

The  resistivity  of  the  small-diameter  wire  used  to  con¬ 
struct  the  coils  is  relatively  high,  increasing  from  33.3 
mfl/m  at  room  temperature  to  40.5  mfl/m  at  75‘^C  (our 
specified  upper  limit  on  temperature  at  the  wire).  The 
additional  length  of  wire  required  to  extend  the  return 
paths  away  from  the  main  cylinders  also  contributed  to 
high  coil  total  resistances.  Since  high  root-mean-squared 
(RMS)  currents  are  required  for  the  desired  imaging  ap¬ 
plications,  the  power  dissipated  by  the  coils  is  significant 
(-1250  W  for  one  coil  @  25A  RMS).  In  a  mock-up  of  one 
quadrant  of  a  coil,  when  a  25 A  RMS  current  was  applied 


with  no  cooling,  the  temperature  of  the  quadrant  rose 
from  10°C  to  100°C  in  2.5  minutes.  With  air  cooling  only, 
a  temperature  rise  from  10°C  to  35°C  in  4  minutes  was 
observed.  A  more  efficient  cooling  system  was  therefore 
required  to  ensure  patient  safety. 

Water  is  one  of  the  most  efficient  and  economical 
cooling  fluids  available.  Its  heat  capacity  (C^  =  4.18 
J/g/°C)  is  approximately  two  times  higher  than  other 
commonly  available  cooling  fluids  such  as  hydrocarbon 
oils  (e.g.,  ethylene  glycol).  Assuming  that  all  heat  gener¬ 
ated  is  transferred  to  the  cooling  water  and  the  maximum 
allowable  temperature  rise  of  the  water  is  5°C,  a  water 
flow  of  3.6  liter/min  is  required  to  remove  1250  W  of 
heat.  Conventionally,  water-cooling  channels  for  gradi¬ 
ent  coils  are  constructed  from  copper  tubing.  For  a  com¬ 
pact  gradient  coil  set,  the  competing  demands  of  space 
limitations  and  the  requirement  of  a  large  cooling  chan¬ 
nel  cross-section  to  ensure  adequate  volume  flow  rates 
necessitated  the  design  of  a  novel  cooling  system.  We 
designed  this  cooling  system  in  separate  modules,  each 
module  corresponding  to  a  quadrant  of  a  single  gradient 
coil.  Each  module  consisted  of  a  skeleton  of  cooling 
channel  spacers  that  were  milled  from  a  3.0-mm-thick 
polycarbonate  sheet.  This  was  covered  on  both  sides 
with  a  very  thin  (0.25-mm-thick)  polycarbonate  film  (GE 
Plastics,  General  Electric  Company,  Pittsfield,  MA)  to 
allow  close  thermal  contact  of  the  wires  with  the  cooling 
water  while  ensuring  electrical  insulation.  The  film  was 
bonded  to  the  skeleton  using  ultraviolet-curing  polyvinyl 
chloride-based  bonding  adhesive  (Loctite  3311,  Loctite 
Corporation,  Mississauga,  Canada).  The  tensile  strength 
(to  break)  of  the  polycarbonate  film  was  rated  at  9000  psi, 
according  to  American  Standards  of  Testing  Materials 
(ASTM)  D882,  and  the  burst  strength  was  200  psi  (ASTM 
D774).  The  modules  were  tested  for  flow  uniformity  and 
resistance  to  leaks  by  flowing  water  at  13  liter/min  for  >5 
hours  before  permanent  incorporation  into  the  gradient 
coil  structure.  After  incorporation  into  the  gradient  coil 
set,  the  modules  were  connected  in  parallel  to  lower  the 
fluid  resistance  of  the  network  of  cooling  channels  and 
thereby  allow  increased  water  flow  rates. 

This  design  allowed  very  close  thermal  contact  be¬ 
tween  the  wires  and  water  in  the  cooling  channels  (while 
maintaining  electrical  isolation)  and  also  provided  a 
large  surface  area  over  which  heat  could  be  transferred  to 
the  water.  The  ability  to  test  the  cooling  channel  integrity 
separate  from  the  gradient  coil  structure  was  vital  in 
ensuring  no  electrical  shorting  due  to  water  leakage.  An 
additional  benefit  of  constructing  the  cooling  system 
from  polycarbonate  is  that  no  electrical  currents  can  be 
created  in  this  material,  and  it  is  therefore  not  susceptible 
to  eddy  current  generation  like  some  copper  systems. 
Since  the  thickness  of  the  cooling  channel  module  was 
only  3.5  mm,  the  maximum  total  wall  thickness  of  the 
gradient  coil  structure  was  less  than  3  cm.  The  assembled 
gradient  coil  structure,  including  RF  coil  and  water  con¬ 
nections,  is  shown  in  Fig.  5. 

RF  Coil 

Several  cylindrical  RF  coil  designs  exist  that  can  be  used 
in  a  transverse  geometry  for  imaging  the  breast  (i.e.,  with 
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FIG.  4.  (a)  Gx  quadrupole  coil  with  216  windings,  (b)  Gy  coil  showing  placement  of  the  “return  wires”  on  a  horizontal  plate. 

Ae  coil  axis  perpendicul^  to  the  main  magnetic  field  proximity  of  the  RF  coil  to  the  gradient  coils  results  in 

f°ij  coils,  when  oriented  vertically,  produce  very  strong  inductive  coupling  and  detuning  of  the  RF 

fields  directed  along  either  the  horizontal  x  direction  coil  (given  that  the  gradient  coil  radius  is  constrained  to 

(e.g.,  saddle  coils  (6))  or  along  the  vertical  direction  (e.g.,  be  as  small  as  possible  for  maximum  efficiency).  To 

solenoidal  coils  (7,  8}).  determine  the  best  RF  coil  design  for  our  gradient  coil 

The  consfiuction  of  an  RF  coil  for  operation  inside  our  set,  two  different  RF  coils  were  constructed  and  tested, 

breast  gradient  coils  is  a  difficult  problem,  in  that  the  The  first  coil  constructed  was  a  quadrature  RF  coil  with 
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FIG.  5.  Assembled  gradient  coll  module  showing  RF  coil  and 
cooling  channels. 

two  separate  coils  providing  the  perpendicular  com¬ 
ponents.  In  this  arrangement,  a  simple  three-loop  sole- 
noidal  coil  provided  the  vertical  mode,  and  a  saddle  coil 
built  onto  the  same  former  provided  the  horizontal  mode. 
The  second  coil  constructed  was  a  16-rung  high-pass 
birdcage  RF  coil  designed  for  operation  in  a  linear  mode 
(recall  that  a  birdcage  coil  cannot  be  used  in  quadrature 
mode  in  the  vertical  orientation  because  one  of  the 
modes  would  be  parallel  to  Bq).  We  chose  this  design 
over  a  saddle-coil  arrangement  because  of  the  superior 
uniformity  of  produced  by  the  birdcage  coil  over  the 
volume  enclosed  by  the  coil.  A  split  cylindrical  RF  shield 
was  constructed  from  fine  copper  screen  and  located 
immediately  inside  the  smallest  diameter  gradient  coil. 
The  coils  were  tested  for  coupling  to  the  gradient  coils. 
The  quality  factors,  Q,  were  measured,  and  the  coil  per¬ 
formances  were  compared. 

Diffusion  Imaging 

Preliminary  diffusion  measurements  were  made  using  a 
short- T2  phantom  consisting  of  NiS04-doped  agarose  gel 
[T2  =  QO  ms  at  1.89  T)  with  a  diffusion-sensitized  spin- 
echo  sequence  (pulsed  gradient  spin-echo  sequence  (9)). 
Apparent  diffusion  coefficient  (ADC)  maps  were  calcu¬ 
lated  from  a  series  of  10  diffusion-weighted  images  with 
diffusion  b  factors  between  0  and  1000  s/mm^.  Each  of 
these  images  was  the  geometrical  average  of  two  images 
with  equal  diffusion-weighting  but  with  opposite  diffu¬ 
sion  gradient  polarity  to  remove  the  effects  of  diffusion¬ 
weighting  cross-terms  from  interaction  between  the  dif¬ 
fusion  and  imaging  gradients  (10).  The  ADC  maps  were 
obtained  under  three  conditions:  (l)  using  the  standard 
whole-body  manufacturer-supplied  gradient  coils  for  all 
imaging  and  diffusion-weighting  gradients;  (2)  using  our 
custom  coils  for  all  imaging  and  diffusion-weighting  gra¬ 
dients,  with  the  gradient  strengths  matched  to  those  of 
the  whole-body  coils  by  applying  an  appropriately  scaled 
current  to  the  coils;  and  (3)  using  our  custom  coils  for  all 
imaging  and  diffusion-weighting  gradients,  with  the  im¬ 
aging  gradient  strengths  matched  to  those  of  the  whole- 
body  coils,  but  taking  advantage  of  the  increased  gradient 
efficiency  to  apply  higher  strength,  shorter  diffusion¬ 
weighting  gradients,  thereby  shortening  the  scan  TE.  In 
the  first  two  experiments,  the  minimum  TE  required  for 


diffusion  weighting  of  b  =  995  s/mm^  was  TE  =  80  ms 
(with  diffusion  gradient  strength  of  22  mT/m).  TR  was  set 
to  600  ms,  the  matrix  size  for  the  scan  was  256  X  128,  and 
the  image  size  was  8X8  cm,  for  a  total  scan  time  of  2.5 
minutes  for  each  image  (5  minutes  for  each  b  factor). 
Diffusion  gradients  were  applied  along  the  readout  axis, 
X,  for  all  experiments.  For  the  third  experiment,  the  TE 
was  decreased  to  TE  =  46  ms  by  increasing  the  strength 
of  the  diffusion-weighting  gradient  pulses  by  a  factor  of 
four  (88  mT/m)  to  obtain  much  larger  b  factors  with 
shorter  diffusion  pulses.  All  other  parameters  (including 
slice  thickness)  were  matched  to  the  first  two  experi¬ 
ments. 

RESULTS 

Coil  Performance 

The  field  produced  by  each  of  the  three  coils  was  calcu¬ 
lated  by  direct  numerical  solution  of  the  Biot-Savart  Law 
for  the  wire  paths,  including  return  wire  contributions. 
Calculated  gradient  efficiency  for  all  three  coils,  together 
with  the  measured  resistance  and  inductance,  are  sum¬ 
marized  in  Table  1.  Contour  plots  of  calculated  gradient 
uniformity  for  the  quadrupole  coil  are  shown  in  Fig.  6. 
Figure  6a  shows  the  percentage  deviation  from  gradient 
uniformity  for  a  y-x  plane  through  the  center  of  the  coils. 
A  similar  plot  is  shown  in  Fig.  6b  for  a  x-z  plane,  trans¬ 
verse  to  the  coil  axis,  at  a  height  equal  to  the  center  of  the 
VOI  (i.e.,  located  5.9  cm  below  the  top  surface  of  the 
coils).  These  plots  show  the  excellent  gradient  unifor¬ 
mity  within  the  VOI  that  is  predicted  for  the  coil.  The 
gradient  strength  decreases  near  the  ends  of  the  coil, 
because  the  finite  length  of  the  coil  causes  a  decrease  in 
field  strength.  Similar  results  were  obtained  for  the 
gradient  coil,  with  appropriate  scaling  for  the  number  of 
wires  and  the  coil  radius.  The  performance  of  the  modi¬ 
fied  target  field  Gy  coil  is  shown  in  Figs.  6c  and  6d.  A 
larger  deviation  from  gradient  uniformity  is  present  for 
this  coil,  compared  to  the  and  G^  coils.  This  problem 
is  inherent  to  any  design  for  a  Gy  coil  in  the  transverse 
orientation,  since  there  is  an  additional  gradient  in  the  y 
direction  due  to  the  natural  field  fall-off  that  occurs  as 
the  ends  of  the  coil  are  approached.  To  create  a  Gy  coil 
with  optimum  gradient  efficiency,  we  would  like  to  max¬ 
imize  \B^\  at  both  ends  of  the  coil,  with  B^  a  maximum  in 
the  -l-z  direction  at  one  end  of  the  coil,  and  a  maximum 
in  the  —  z  direction  at  the  opposite  coil  end.  The  field 
fall-off  inherent  in  this  coil  orientation  works  against  the 
gradient  we  desire  to  create,  since  it  necessarily  results  in 
a  decrease  in  the  magnitude  of  B^  as  the  coil  ends  are 
approached.  Instead  of  working  against  the  natural  field 
fall-off  toward  the  ends  of  the  coil,  we  have  designed  a 


Table  1 

Coil  Performance  Characteristics 


Mean 

diameter 

(cm) 

Efficiency 

(mT/m/A) 

Inductance 

(ixH) 

Resistance 

(n) 

(quadrupole) 

21.1 

2.9 

754 

2.3 

Gy  (fingerprint) 

19.5 

1.7 

829 

2.6 

(quadrupole) 

18.1 

3.3 

1024 

1.9 
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FIG.  6.  Contour  plots  showing  percentage  deviations  from  gradient  uniformity  for  the  coil  over  (a)  an  axial  plane  (x-y)  showing  the 
entire  VOI  with  the  top  edge  of  the  plot  corresponding  to  the  top  surface  of  the  coils.  The  outer  circumference  of  the  VOI  is  at  6.75  cm, 
as  indicated  by  dashed  lines,  (b)  A  coronal  plane  (x-z)  located  at  the  center  of  the  defined  VOI,  i.e.,  at  ~3  cm  above  the  geometric  coil 
center.  In  Fig.  6a,  the  contours  are  symmetric  around  y  =  0  cm  with  some  small  deviations  due  to  return-wire  contributions,  (c,  d)  The 
corresponding  plots  for  the  constructed  Gy  coil  shown  in  Fig.  4b.  (e,  f)  The  corresponding  plots  for  a  conventional  transverse  coil  with 
return  wires  on  the  cylinder,  as  in  Fig.  3b. 


coil  to  produce  0  at  the  upper  end  of  the  coil,  and  a 
maximum  in  toward  the  bottom  end,  similar  to  the 
asymmetric  head  coil  design  by  Roemer  (11).  We  used 
the  target  field  approach  to  optimize  the  rate  at  which  the 
field  falls  off  toward  the  upper  end  of  the  coil  to  produce 
a  linear  variation  in  B^.  Although  the  region  of  acceptable 
gradient  uniformity  for  our  truncated  target  field  coil  is 
not  as  large  as  for  the  quadrupole  coils,  it  compares 
favorably  to  the  performance  of  the  full  target  field  de¬ 
sign.  We  show  in  Figs.  6e  and  6f  the  corresponding  plots 
from  a  complete  40  loop  per  side  target  field  coil  with  the 
same  dimensions  as  our  modified  coil  (i.e.,  with  return 
wire  paths  at  the  bottom  of  the  cylinder  as  in  Fig.  3b).  The 
region  of  acceptable  gradient  uniformity  for  our  trun¬ 
cated  target  field  coil  is  much  larger  than  for  the  com¬ 
plete  target  field  coil.  In  addition,  the  gradient  efficiency 
is  better  for  our  truncated  coil  by  a  factor  of  1.7,  due  to 


the  larger  number  of  windings  that  could  be  placed  on 
the  cylinder  for  our  coil. 

To  verify  the  field  simulations,  gradient  uniformity 
was  tested  by  acquiring  images  of  grid  phantoms.  Devi¬ 
ations  from  gradient  uniformity  are  manifest  in  the  im¬ 
ages  as  geometric  distortions.  The  phantoms  consisted  of 
acrylic  cylinders  filled  with  square  plastic  gridding  and 
CuS04-doped  0.9%  saline  solution.  The  grid  spacing  was 
1.7  X  1.7  cm,  whereas  the  phantoms  had  an  outer  diam¬ 
eter  of  12.5  cm  and  a  height  of  15  cm  (extending  3  cm 
higher  than  the  top  surface  when  placed  inside  the  coils). 
The  coil  performance  was  tested  on  a  clinical  GE  Signa 
imager,  using  our  coils  for  all  imaging  gradients  (with 
gradient  strengths  matched  to  the  whole-body  coils)  in  a 
standard  spin-echo  pulse  sequence.  In  general,  local  gra¬ 
dient  coils  often  require  modification  of  the  compensa¬ 
tion  network  for  the  gradient  coil  amplifiers  since  the 


398 


Maier  et  al. 


coil  inductances  are  quite  different  from  the  whole-body 
coils.  Because  the  inductances  of  our  breast  coils  were 
very  similar  to  the  whole-body  coils,  no  modification  of 
the  network  was  necessary.  Images  of  a  grid  phantom  are 
shown  in  Fig.  7  for  planes  corresponding  approximately 
to  the  theoretical  field  plots  shown  in  Figs.  6a  and  6b  and 
are  evidence  of  the  good  gradient  uniformity  achieved 
over  most  of  the  VOL  These  images  were  acquired  using 
a  spin-echo  sequence  with  TE  =  20  ms,  TR  =  500  ms, 
matrix  size  =  256  X  256,  and  2  NEX.  The  image  size  is 
16  X  16  cm,  and  the  VOI  vertical  extent  is  marked  with 
arrows. 

The  first  RF  coil  we  constructed  was  a  quadrature  coil 
consisting  of  a  solenoidal  coil  to  provide  the  vertical 
mode  and  a  saddle  coil  to  provide  the  horizontal  coil. 
Unfortunately,  the  vertical  mode  was  inductively  cou¬ 
pled  to  both  the  y-gradient  coil  and  to  the  RF  shield  and 
could  not  be  tuned  to  the  proton  resonant  frequency  at 
1.5  T.  We  therefore  chose  to  construct  a  coil  that  would 
operate  in  linear  mode,  with  only  a  horizontal  com¬ 
ponent  present.  We  constructed  a  birdcage  coil  that  cre¬ 
ates  a  horizontal  mode  field  with  a  high  degree  of 
homogeneity  over  a  large  volume.  The  unloaded  Q  of  this 
coil  (when  located  inside  the  gradient  coil  structure)  was 
80,  whereas  the  Q  for  loading  with  a  0.9%  saline  phan¬ 
tom  was  29. 

Although  the  gradient  strengths  required  for  conven¬ 
tional  imaging  (10  mT/m)  can  be  generated  at  less  than  6 
A  using  these  coils,  higher  currents  are  required  for  more 
demanding  applications  such  as  flow  imaging  and  diffu¬ 
sion  imaging.  To  test  the  effectiveness  of  the  water  cool¬ 
ing  system,  a  DC  current  was  applied  to  the  coil  with  the 
highest  total  resistance,  Gy,  and  the  temperature  rise  was 
measured  for  different  water  flow  rates  by  monitoring  a 
thermocouple  located  on  the  surface  of  the  wires.  The 
temperature  was  monitored  until  equilibrium  was 


achieved.  Measured  temperature  rises  are  shown  in  Fig.  8 
for  DC  currents  between  0  and  22  A  and  a  cooling  water 
flow  rate  of  8  liter/min.  Similar  results  were  obtained  for 
other  flow  rates,  at  which  the  equilibrium  temperature 
increased  with  decreasing  flow  rates  but  did  not  change 
significantly  from  the  results  shown  in  Fig.  8  for  higher 
flow  rates.  This  behavior  has  also  been  demonstrated  in 
purely  cylindrical  gradient  coils  and  was  modeled  by 
Chu  and  Rutt  (12)  using  an  analysis  based  on  fundamen¬ 
tal  heat  transfer  theory. 

Diffusion  Imaging 

ADC  maps  for  a  T2  =  60  ms  NiS04-doped  agar  phantom 
are  shown  in  Fig.  9;  Fig.  9a  shows  the  results  obtained 
using  the  GE  gradient  coils,  and  Fig.  9b  shows  the  results 
obtained  using  our  gradient  coils  with  diffusion  gradient 
strengths  four  times  higher  than  that  of  the  GE  coils.  The 
mean  ADC  for  all  three  experiments  was  the  same  within 
experimental  error:  (1.8  ±  0.3)  X  10“^  mm^/s  for  the  GE 
coils  (±30%  peak-to-peak),  (1.8  ±  0.3)  X  10"^  mm^/s  for 
our  gradient  coils  with  diffusion  gradient  strengths 
matched  to  the  GE  coils,  (1.9  ±  0.2)  X  10"^  mm^/s  for  our 
gradient  coils  with  diffusion  gradient  strengths  four 
times  higher  than  that  of  the  GE  coils  and  TE  =  46  ms.  As 
expected,  the  observed  signal-to-noise  ratio  was  in¬ 
creased  by  nearly  a  factor  of  2  for  the  TE  =  46  ms  images 
compared  to  the  corresponding  h-factor  images  at  TE  = 
80  ms. 

DISCUSSION 

The  performance  and  homogeneity  of  the  gradient  coils 
were  in  excellent  agreement  with  theoretical  predictions. 
The  volume  of  acceptable  gradient  uniformity  (<20% 
deviation)  extended  over  most  of  our  defined  VOI  for  all 


a  b 

FIG.  7.  Spin-echo  images  of  a  grid  phantom  acquired  (a)  for  an  axial  plane  (x-y)  and  (b)  for  a  coronal  plane  (x-z)  with  the  specialized  breast 
quadrupole  gradient  coils.  The  input  current  strength  required  for  the  least  efficient  (G^)  coil  was  approximately  14  times  less  than  that 
required  for  the  corresponding  GE  coil. 
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FIG.  8.  Temperature  response  of  assembled  gradient  colls  to  a 
DC  current  applied  to  the  coil  with  the  highest  resistance  (Gy) 
monitored  at  the  surface  of  the  Gy  coil  where  the  winding  density 
was  greatest. 


three  gradient  coils.  For  this  region,  the  theoretical  plots 
of  can  be  used  to  generate  a  gradient-warping  algo¬ 
rithm  for  the  coils  (13)  to  correct  spatial  distortions  in  the 
images.  Additionally,  the  known  deviations  from  gradi¬ 
ent  uniformity  can  be  used  to  correct  calculations  of 
diffusion-weighting  b  factors  (which  will  vary  spatially 
across  the  VOI  for  our  coils).  Unfortunately,  greater  de¬ 
viations  from  uniformity  could  not  be  avoided  for  the  top 
1  cm  of  the  VOI  near  the  upper  edges  of  the  coils  for  our 
coil  designs.  Deviations  from  uniformity  of  more  than 
20%  lead  to  large  image  distortions  and  likely  cannot  be 
satisfactorily  corrected  using  conventional  schemes. 
Many  breast  lesions  occur  close  to  the  chest  wall;  to  be 
clinically  useful,  it  is  important  to  image  accurately  as 
close  to  the  chest  wall  as  possible.  Several  new  gradient 
coils  with  improved  gradient  uniformity  at  the  top  edge 
of  the  coil  have  been  reported  recently  (14,  15).  Future 
improvements  on  our  coil  design  will  focus  on  extending 
the  region  of  usable  gradient  beyond  the  top  edge  of  the 
coil  (16)  as  far  as  possible.  For  our  preliminary  clinical 
work  with  these  coils,  we  will  limit  our  study  to  patients 
with  known  lesions  occurring  in  a  region  of  the  breast 
that  falls  within  the  20%  uniformity  region  of  our  gradi¬ 
ent  coils.  This  approach  will  allow  us  to  implement 
pulse  sequences  under  ideal  conditions,  (a)  to  determine 
the  best  parameters  for  visualizing  breast  lesions,  (b)  to 
determine  the  minimum  required  gradient  strengths  and 
uniformities,  and  (c)  to  make  an  initial  assessment  of  the 
clinical  utility  of  this  approach. 

In  future  coil  designs,  it  may  be  possible  to  extend  the 
region  of  usable  gradient  farther  into  the  cylinder  in  the 
~y  direction  by  recessing  the  coils  in  the  patient  table  or 
by  constructing  a  dedicated  breast  imaging  table,  thereby 
allowing  higher  aspect  ratio  coils.  In  this  design,  the  coil 
assembly  would  be  constructed  to  allow  positioning  on 


either  side  of  the  table,  depending  on  which  breast  was 
being  imaged.  A  fiberglass  or  Kevlar  sheath  covering  the 
coil  assembly  and  running  the  length  of  the  patient  table 
would  serve  the  dual  functions  of  isolating  the  patient 
from  the  electrical  and  cooling  circuits  and  providing  an 
additional  form  of  mechanical  fixation  to  secure  the  coils 
on  the  table.  Because  the  coils  are  not  torque-balanced,  it 
will  be  necessary  to  anchor  the  structure  firmly,  both 
from  the  point  of  view  of  vibration-induced  imaging  ar¬ 
tifacts  and  of  patient  safety. 

The  high  gradient  efficiencies  of  these  coils  will  allow 
much  higher  spatial-resolution  imaging  than  can  cur¬ 
rently  be  achieved  in  a  clinical  setting,  and  additionally, 
will  provide  very  large  amplitude  motion-sensitizing  gra¬ 
dients,  thus  facilitating  diffusion  imaging  (or  flow  imag¬ 
ing)  in  small  blood  vessels  with  short  echo  times.  T2  has 
been  measured  for  human  breast  tissue,  with  reported 
values  ranging  between  60  and  150  ms  (17, 18).  Recently, 
we  reported  T2  approximately  70  ms  for  human  breast 
tumor  (MCF7  human  breast  carcinoma)  growing  in  a 
mouse  mammary  fat  pad  (19).  More  fibrous  tumors  can 
have  significantly  shorter  T2S.  For  this  reason,  initial 
signal  loss  at  h  =  0  s/mm^  due  to  T2  weighting  of  the 
image  can  severely  limit  the  range  of  b  factors  for  which 
images  can  be  obtained,  which  has  the  result  of  limiting 
the  precision  of  the  ADC  measurement.  Our  approach  is 
to  shorten  TE  while  maintaining  diffusion  weighting  by 
increasing  diffusion  gradient  amplitudes.  This  gives  a 
substantial  increase  in  initial  signal  intensity  for  h  =  0 
s/mm^,  allowing  high-quality  images  to  be  obtained  for 
much  higher  b  factors.  Even  in  cases  in  which  signal  loss 
from  T2  weighting  is  not  so  severe  as  to  limit  ADC  pre¬ 
cision,  the  ability  to  obtain  high  diffusion  weighting  with 
short  TE  reduces  dramatically  the  T2  contrast  in  an  im¬ 
age.  This  will  allow  acquisition  of  a  diffusion-weighted 
image  in  which  the  contrast  is  determined  mainly  by 
tissue  ADC  rather  than  by  a  mix  of  ADC  and  T2. 

A  secondary  important  benefit  that  can  be  gained  using 
high-efficiency  gradient  coils  is  that  the  ramp  time  re¬ 
quired  to  attain  imaging  gradient  strengths  (10-25 
mT/m)  is  decreased  from  that  required  by  whole-body 
coils  by  a  factor  of  the  efficiency  (since  the  current  slew 
rate  is  similar  to  whole-body  coils,  and  the  required  peak 
current  is  decreased).  This  will  enable  echo-planar  imag¬ 
ing  (and  other  fast  imaging  techniques)  using  conven¬ 
tional  gradient  amplifiers.  Echo-planar  imaging  tech¬ 
niques  have  been  proven  to  be  much  less  susceptible  to 
motion  artifacts  than  conventional  imaging  sequences 
(20)  and  have  recently  been  applied  to  diffusion  imaging 
in  vivo  (21)  with  good  results. 

For  rapid  imaging  pulse  sequences,  and  for  high-am- 
plitude  bipolar  pulses,  very  rapid  changes  in  gradient 
polarity  are  required.  Although  the  magnetic  field  ampli¬ 
tudes  drop  off  quickly  in  the  region  above  the  coils, 
magnetic  fields  are  also  generated  by  the  return  wires  on 
the  horizontal  plate,  and  excitation  thresholds  for  cardiac 
stimulation  must  be  considered.  As  an  upper  limit  for 
one  coil,  \B\  was  calculated  for  the  coil  at  various 
distances  above  the  return  wires  at  the  location  of  the 
highest  wire  density.  The  maximum  dl/dt  from  our  gra¬ 
dient  coil  amplifiers  is  200  A  in  184  /uis.  This  produces 
dB/dt  =  156T/satl  cm  above  the  coil  return  wire  plate 
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FIG.  9.  ADC  maps  of  a  Tg  =  60  ms  agar  phantom  obtained  (a)  using  the  whole-body  coils  to  provide  ali  diffusion  and  imaging  gradients, 
TE  =  80  ms,  and  (b)  using  the  breast  gradient  coiis  to  provide  all  gradients,  where  the  diffusion  gradient  strengths  were  four  times  the 
whole-body  gradient  strengths  and  TE  =  46  ms.  In  all  cases,  the  ADC  maps  were  constructed  from  10  images  with  b  =  0  to  b  =  1000 
s/mm^. 


and  dB/dt  =  50  T/s  at  10  cm,  compared  to  a  recom¬ 
mended  upper  limit  of  400  T/s  for  cardiac  stimulation 
(22).  The  magnitude  of  these  fields  is  relatively  small  at 
regions  of  concern  with  respect  to  magnetic  stimulation 
such  as  the  heart;  however,  peripheral  nerve  stimulation 
has  heen  reported  at  dB/dt  as  low  as  40  T/s  in  some 
patients  (22).  Theoretical  field  calculations  will  allow 
prediction  of  \B\  and  dB/dt  for  all  pulse  sequences  before 
implementation,  and  pulse  sequences  will  be  tailored  to 
adhere  to  maximum  dB/dt  and  peak  \B\  guidelines.  Ad¬ 
ditionally,  if  peripheral  nerve  stimulation  is  reported  for 
individual  patients,  dB/dt  will  be  reduced  below  the 
stimulation  threshold  by  increasing  gradient  switching 
times. 


SUMMARY 

We  have  constructed  a  three-axis  gradient  coil  set  for 
imaging  human  breast  tumors.  These  coils  create  gradi¬ 
ent  efficiencies  of  up  to  3.3  mT/m/A  with  relatively  low 
inductance  (<1000  /xH),  and  good  gradient  uniformity 
over  most  of  the  volume  required  for  breast  imaging. 
Preliminary  diffusion  results  were  obtained  in  a  water 
phantom  for  b  factors  as  high  as  ~1000  s/mm^  for  TE  = 
46  ms.  A  compact,  efficient  water-cooling  system  was 
presented,  which  maintained  coil  temperatures  at  or  be¬ 
low  room  temperature,  even  for  RMS  currents  as  high  as 
20  A.  This  gradient/cooling  module  will  allow,  for  the 
first  time,  clinical  breast  imaging  using  diffusion  and 
microvascular  flow-weighted  protocols  with  high  spatial 
resolution. 
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,  „  ,  MAGNETIC  RESONANCE  IN  MEDICINE 

Constrained  Length  Minimum  Inductance  Gradient 
Coil  Design 

Blaine  A.  Chronik,  Brian  K.  Rutt 


A  gradient  coil  design  algorithm  capable  of  controlling  the 
position  of  the  homogeneous  region  of  interest  (ROI)  with 
respect  to  the  current-carrying  wires  is  required  for  many 
advanced  imaging  and  spectroscopy  applications.  A  modified 
minimum  inductance  target  field  method  that  allows  the 
placement  of  a  set  of  constraints  on  the  final  current  density 
is  presented.  This  constrained  current  minimum  inductance 
method  is  derived  in  the  context  of  previous  target  field  meth¬ 
ods.  Complete  details  are  shown  and  all  equations  required 
for  implementation  of  the  algorithm  are  given.  The  method  has 
been  implemented  on  computer  and  applied  to  the  design  of 
both  a  1:1  aspect  ratio  (lengthidiameter)  central  ROI  and  a  2:1 
aspect  ratio  edge  ROI  gradient  coil.  The  1:1  design  demon¬ 
strates  that  a  general  analytic  method  can  be  used  to  easily 
obtain  very  short  gradient  coil  designs  for  use  with  specialized 
magnet  systems.  The  edge  gradient  design  demonstrates  that 
designs  that  allow  imaging  of  the  neck  region  with  a  head¬ 
sized  gradient  coil  can  be  obtained,  as  well  as  other  applica¬ 
tions  requiring  edge-of-cylinder  regions  of  uniformity. 

Key  words:  gradient  coils;  target  field;  current  constraints; 
minimum  inductance;  constrained  length. 

INTRODUCTION 

Advanced  MR  imaging  and  spectroscopy  methods  often 
require  high  performance  gradient  hardware,  where  per¬ 
formance  is  typically  characterized  by  a  combination  of 
gradient  efficiency  (gradient  strength  per  unit  current) 
and  switching  speed.  The  design  of  gradient  coils  cus¬ 
tomized  to  a  particular  application  is  a  method  of  realiz¬ 
ing  gains  in  both.  One  defining  characteristic  of  custom¬ 
ized  gradient  coils  is  that  they  conform  closely  to  the 
anatomical  region  of  interest  and  although  this  allows 
increased  performance,  subject  access  to  the  usable  re¬ 
gion  of  the  coil  becomes  an  important  and  nontrivial 
design  consideration. 

Gradient  coil  design  methods  can  be  divided  into  two 
general  categories:  “coil  space”  versus  “reciprocal  space” 
techniques.  Coil  space  methods  iteratively  adjust  the  po¬ 
sitions  of  wires  on  the  gradient  coil  former  until  the 
design  conforms  to  some  desired  performance  specifica¬ 
tions.  The  optimal  wire  configuration  is  chosen,  using  a 
numerical  optimization  technique,  from  the  multitude  of 
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wire  configurations  that  span  the  domain  of  allowable 
patterns.  The  method  of  Conjugate  Gradient  Descent  (1, 
2)  and  the  method  of  Simulated  Annealing  (3-6)  are 
examples  of  different  optimization  techniques  applied  to 
this  problem.  Although  these  methods  in  principle  allow 
complete  freedom  in  the  wire  pattern  and  can  be  applied 
to  any  geometry,  they  are  in  practice  limited  due  to  the 
computational  intensity  of  the  methods.  They  have  been 
most  effectively  applied  to  the  problem  of  positioning 
coaxial  circular  wire  loops  in  the  design  of  longitudinal 
(GJ  gradient  coils;  furthermore,  when  applied  to  the 
design  of  cylindrical  transverse  gradient  coils  (G^  or  G^), 
the  functional  form  of  the  wire  shapes  have  typically 
been  predetermined  to  depend  on  only  a  few  parameters 
so  that  the  number  of  degrees  of  freedom  remains  prac¬ 
tical.  Advantages  of  these  methods  include  the  absence 
of  an  approximation  stage  in  sampling  a  continuous  cur¬ 
rent  density  with  a  finite  number  of  wires,  and  most 
importantly  the  ease  with  which  they  incorporate  direct 
constraints  on  the  final  wire  pattern,  thereby  allowing 
the  design  of  restricted  length  gradient  coils.  The  first 
advantage  is  primarily  of  concern  only  when  the  current 
density  is  to  be  sampled  with  a  small  number  of  wires, 
since  field  variations  due  to  finite  wire  spacing  become 
insignificant  at  distances  larger  than  that  spacing. 

Reciprocal  space  methods  are  characterized  by  the  ex¬ 
pansion  of  all  quantities  of  interest  in  terms  of  orthogonal 
functions  and  solving  for  the  coefficients  of  the  corre¬ 
sponding  current  density  expansion.  The  current  density 
is  then  approximated  by  a  finite  set  of  wires,  typically 
using  stream  function  methods  (7).  Because  of  the  diffi¬ 
culty  in  obtaining  orthogonal  function  expansions  for 
general  geometries,  in  practice  these  techniques  are  lim¬ 
ited  to  simple  geometries  (planes  (8,  9),  cylinders  (10, 
11),  spheres  (12),  ellipses  (13),  hyperbolas  (14));  further¬ 
more,  they  often  lack  control  over  the  final  current  den¬ 
sity  that  is  attainable  with  coil  space  methods  (15).  Once 
quantities  of  interest  have  been  expanded  and  the  coef¬ 
ficients  of  the  current  density  expansion  formally  solved 
for,  computer  implementation  of  these  methods  is  simple 
and  current  densities  can  be  obtained  very  quickly.  Re¬ 
ciprocal  space  methods  can  yield  optimal  designs  by 
using  constrained  Lagrange  multiplier  or  least  squared 
error  methods  to  solve  for  the  expansion  coefficients  (16). 

In  this  paper  we  address  the  problem  of  including 
direct  control  over  current  density  in  a  reciprocal  space 
design  method.  We  present  the  details  of  a  modification 
of  Turner’s  minimum  inductance  target  field  method  (17) 
that  allows  explicit  constraint  of  the  extent  of  the  final 
current  density  as  well  as  control  over  the  position  of  the 
uniform  gradient  with  respect  to  the  current  density 
while  minimizing  coil  inductance.  These  modifications 
yield  a  design  method  with  the  generality  claimed  for 
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coil  space  methods  and  the  efficiency  of  reciprocal  space 
techniques. 


THEORY 

In  any  magnetic  design  problem  the  most  important  step 
is  the  establishment  of  a  relationship  between  the  current 
density  function  and  the  desired  magnetic  field  compo¬ 
nents.  For  gradient  coils  used  in  MRI,  we  are  only  con¬ 
cerned  with  the  z-component  of  the  magnetic  field,  B^{r, 
<j>,  z];  furthermore,  if  we  are  considering  axially  aligned 
cylindrical  gradient  coils  (define  the  coil  axis  to  lie  along 
the  z  direction),  only  the  azimuthal  component  of  current 
density,  J^[^,  z),  contributes  to  BJ^r,  (/>,  z)  and  therefore 
the  problem  reduces  to  obtaining  the  connection  between 
these  two  quantities.  There  are  several  ways  to  express 
this  relationship  but  one  that  has  proved  especially  fruit¬ 
ful  is  the  “target  held’'  formalism  first  described  by 
Turner  (10).  This  formalism  is  characterized  by  the  ex¬ 
pansion  of  B^[r,  (f),  z)  in  cylindrical  harmonics  (18): 


and 

(+77  l+eo 

B^(c,  k)=  d<i>  e-™'''  dz  4>,  z)  [5] 

~7T  *1-00 

with  c  the  radius  of  a  cylindrical  surface  (c  <  a),  called 
the  “target  cylinder.”  It  is  useful  to  think  of  B^(c,  k)  as  an 
apodised  version  of  j^{k),  and  o/"[c,  k)  acting  as  the 
apodising  function  with  limiting  forms; 

k)  ~  k^oo  [6a] 

k  =  0)  =  —  [6b] 
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Knowing  that  B^{c,  k)  and  j^[k)  are  related  via  a  simple 
product  translates  into  a  convolution  relationship  be¬ 
tween  B^[Cy  </),  z)  and  z): 

B^{c,  4>,  z)  =  0,(c,  z)®J^{4),  z)  [7] 
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where  and  are  modified  Bessel  functions,  a  is  the 
radius  of  the  cylinder  on  which  the  current  is  con- 
strained  to  flow,  and 


is  the  2D  Fourier  Transform  of  z)  in  cylindrical 
coordinates.  All  target  field  based  methods  express  quan¬ 
tities  of  interest  as  functions  of  j^[k)  and  then  directly 
solve  for  the  j^[k]  that  yields  the  desired  gradient  coil. 
Methods  only  differ  in  terms  of  the  physical  quantities 
explicitly  included  in  the  design  process  and  how  the 
required  j^[k)  is  solved  for.  It  should  be  emphasized  that 
all  target  field  methods  are  direct  inverse  solutions  to  the 
gradient  coil  design  problem  and  that  once  j^[k]  is 
known,  the  current  density  obtained  via  inverse  transfor¬ 
mation  of  j^{k)  will  give  the  desired  gradient  coil,  i.e., 
one  which  produces  exactly  the  specified  field. 

“Original”  Target  Field 

The  original  target  field  method  developed  by  Turner 
(10)  makes  use  of  the  fact  that  the  2D  Fourier  transform  of 
B^{r,  (j),  z)  is  very  simply  related  to  j^{k]: 

BTic,  k)  =  k)j^{k)  [3] 

where 

a)-(c,  k)  =  i-  tioCt)\k\U\27rkc\)KU\27Tka\)  [4] 


where  n(c,  z)  is  the  inverse  transform  of  co"”(c,  k).  D(c,  z) 
is  approximately  Lorentzian  in  shape  with  full  width  at 
half-maximum  proportional  to  (a-c);  this  means  that  the 
magnetic  field  has  the  same  functional  form  as  the  azi¬ 
muthal  current  density  but  “smeared”  in  the  z  direction 
by  a  Lorentzian-like  function  that  gets  broader  as  the 
target  cylinder  moves  further  away  from  the  current- 
carrying  cylinder.  In  other  words,  the  longitudinal  com¬ 
ponent  of  magnetic  field  is  a  “low-pass  filtered”  version 
of  the  azimuthal  component  of  current  density. 

To  obtain  a  gradient  coil,  the  desired  field  on  the  target 
cylinder  must  be  specified  (i.e.,  the  “target  field”).  Gen¬ 
erally,  we  write; 

B^{c,  (/),  z)  =  g(z)/(c,  (jy)  [8] 

where  g(z)  is  the  field  modulation  along  the  z  direction 
and/(c,  (f))  describes  the  radial  and  azimuthal  variation  of 
B^.  For  longitudinal  (GJ  gradient  coils  we  would  have 
y(c,  <^)  =  constant  and  the  regions  of  uniform  gradient 
correspond  to  regions  over  which  the  first  derivative  of 
g(z)  is  approximately  constant.  For  transverse  gradient 
coils  (G^  or  Gy)  we  would  have: 


f(c,  </))  =  cGx  cos  (j) 

[GJ 

[9a] 

f(Cy  (^)  =  cGy  sin  (f) 

[Gy] 

[9b] 

and  regions  of  uniform  gradient  correspond  to  regions 
over  which  g(z)  is  approximately  constant. 

Taking  a  G^  coil  as  an  example,  we  combine  Eqs.  [8] 
and  [9a],  substitute  the  resulting  expression  for  B^{c,  k) 
into  Eq.  [3]  and  solve  for  j^ik)  to  get: 


■in\  cG^7r-S{g(z)} 

0)  (c,  k) 

[10a] 

i^{k)  =  0;  for  m  1 

[10b] 

where  ^{g(z)}  is  the  ID  Fourier  transform  of  g(z).  Once 
g(z)  has  been  specified  over  all  z  points,  the  ^(k)  func- 
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tion  is  easily  arrived  at  using  Eq.  [10a]  and  z]  ob¬ 
tained  by  inverse  transformation.  The  z)  current 
density  component  is  obtained  from  z]  by  imposing 
continuity  of  J{<^,  z)  (7,  19).  It  should  be  noted  that  al¬ 
though  the  forms  in  this  paper  are  general  in  that  they 
include  the  complete  range  of  possible  m  values,  the 
expressions  applied  to  longitudinal  and  transverse  gra¬ 
dient  coils  greatly  simplify  to  include  only  the  m  =  0  and 
m  =  ±1  contributions  respectively. 

Minimum  Inductance  Target  Field 

The  minimum  inductance  target  field  method  represents 
a  rather  different  and  more  po’werful  way  of  solving  for 
f^ik)  in  that  we  now  determine  the  j^[k]  that  minimizes 
some  property  of  the  gradient  coil  (typically  inductance 
or  power  dissipation)  subject  to  the  constraint  that  be 
equal  to  specified  values  [B^^)  at  a  finite  set  of  coordi¬ 
nates  (r^,  <f)^,  z^).  The  result  of  the  analysis  is  still  j^{k) 
(the  transform  of  which  gives  the  required  current  den¬ 
sity)  and  Eq.  [3]  still  holds;  however,  j^[k]  is  solved  for 
differently. 

To  set  a  context  for  the  extension  described  in  this 
paper,  we  review  the  form  of  the  minimum  inductance 
method,  first  described  by  Turner  (17),  It  is  an  applica¬ 
tion  of  the  method  of  constrained  optimization  using 
Lagrange  multipliers  (19).  A  functional,  U{j^{k)},  is 
formed  that  includes  the  function,  Flj^{k]]  to  be  opti¬ 
mized  (typically  F  is  either  inductance,  L,  or  power  dis¬ 
sipation,  P,  both  of  which  are  expressible  as  functions  of 
the  current  density)  and  the  field  constraint  equations 
multiplied  by  constants,  the  Lagrange  multipliers: 

rnf^ik)}  =  F{j^(k)}  +  2;  ,  </.„ ,  z„)  -  B J  [11] 

n 

In  a  manner  analogous  to  the  optimization  of  a  func¬ 
tion  of  a  single  variable,  the  derivative  of  U  with  respect 
to  j^[k)  is  taken  and  set  equal  to  zero.  We  then  solve  that 
equation  for  /J(k).  For  the  specific  case  of  minimizing 
coil  inductance,  we  first  write  L  as  a  function  of  j^[k) 
(17,  20): 


[12] 


r^{\2Trka\)K',„i\2TTka\)j^{ky 


where  I  is  the  input  current  required  to  produce  the 
surface  current  density,  then  note  the  form  of  the  field 
constraint  equations: 


<l>„,  zj  -B^„  =  0  [13a] 


+  00 

m=:~00  •'  -00 

[13b] 

IJl2irkrJ)K’„(l2iTkal)j^(k)  -B^„  =  0 


Differentiating  U  composed  of  Eqs.  [12]  and  [13b]  and 
solving  for  j^[k)  yields: 

JV 

mk)  =  E  A„a”(lc)  [14] 

71  =  1 


I^i\2Trkr„\) 

rj\27rka\) 


[15] 


Once  the  set  of  is  known,  Eqs.  [14]  and  [15]  give  the 
required  j^{k).  The  A^’s  are  determined  by  substituting 
j^{k)  as  expressed  by  Eq.  [14]  into  the  N  constraint  equa¬ 
tions  of  Eq.  [13b].  This  can  be  Avritten  as  a  matrix  equa¬ 
tion  for  the  A^’s: 


[Mii][An]  —  [Bzn]  [16] 

where  the  elements  of  matrix  are  known  functions  of 
the  constraint  coordinates  and  are  given  as  the  first  row 
of  Table  1.  The  design  algorithm  is  then  simply  to  choose 
the  N  field  constraints,  calculate  the  elements  of 
from  the  constraint  coordinates  and  the  equations  shovra 
in  Table  1,  solve  Eq.  [16]  for  the  set  of  {A^},  compose  j^[k] 
via  Eq.  [14]  and  inverse  transform  to  get  /^((/),  z). 


Constrained  Current  Minimum  Inductance  Target  Field 

To  add  explicit  control  over  the  length  (and  boundaries) 
of  the  final  current  density  to  the  design  algorithm,  we 
have  modified  Turner’s  minimum  inductance  method  to 
include  a  set  of  current  constraints.  This  type  of  control 
over  the  design  is  very  important  when  the  region  of  good 
uniformity  is  required  to  be  near  the  extremes  of  the 
physical  coil,  as  well  as  when  the  total  gradient  coil 
length  needs  to  be  minimized. 

The  first  form  of  current  constraint  that  one  would 
want  to  impose  upon  a  design  is  a  “closure”  constraint 
that  prevents  any  current  density  from  crossing  the 
boundaries  of  a  specified  area  on  the  surface  of  the  cy¬ 
lindrical  former.  This  constraint  manifests  itself  in  the 
final  coil  design  as  forcing  all  conducting  loops  within 
the  specified  region  to  be  closed.  A  simple  form  of  this 
constraint  would  be  to  require  that  no  cmrent  density 
cross  the  azimuthal  lines  defined  by  ±Zg.  Analytically 
this  is  equivalent  to  requiring  that  the  integral  of  /^(4>,  z) 
on  the  line  along  any  (p  =  <l>q  from  -z^  to  +Zg  be  zero: 

r  +Z7 

Zq)  =  I  IMq’  =  0  [17] 

•'-0, 

For  longitudinal  gradient  coils  this  constraint  does  not 
need  to  be  applied  because  the  z  component  of  current 
density  will  always  be  zero.  For  transverse  gradient  coils 
the  azimuthal  variation  of  /^(4>,  z)  is  always  sinusoidal; 
therefore,  the  integral  constraint  only  needs  to  be  applied 
over  a  single  line  from  -z^  to  -l-Zg  (i.e.,  along  constant 
</>  =  <#>,)• 

The  second  type  of  current  constraints  necessary  to 
force  the  current  density  to  remain  contained  within 
some  region  of  the  coil  surface  are  direct  constraints  on 
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Table  1 

Analytic  Forms  for  the  Calculation  of  Matrix  Elements  and  Subsequent  Determination  of  Lagrange  Multipliers  in  the  Constrained 
Current  Minimum  Inductance  Design  Method 


M,,{n 


M^2{n,  p 


MiaCn,  g') 


jn  =  -oo  oo 

f»l  =  ~oo 
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+  00  M 

12 
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w^(k) 
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the  values  of  the  azimuthal  current  density  over  a  set  of 
coordinates: 


^p)  J^p’f  p  —  1:P  [18] 

Azimuthal  current  constraints  of  this  form  have  been  intro¬ 
duced  previously  (21);  however,  it  should  be  noted  that  only 
the  use  of  both  the  above  forms  of  current  constraint  (closure 
and  point  azimuthal  constraints)  provides  a  sufficient  and 
necessary  condition  for  current  density  to  be  zero  beyond  ±Zg. 
The  closure  constraint  alone  does  not  disallow  closed  current 
loops  outside  the  bounded  region  whereas  constraining  the 
azimuthal  current  density  to  be  zero  does  not  prohibit  non¬ 
zero  longitudinal  current  density,  z)  at  the  P  constraiut 
points;  this  will  be  made  clear  by  example  in  the  Results 
section. 

Mathematically  the  current  constraints  represent  a  sim¬ 
ple  extension  of  the  minimum  inductance  method  outlined 
above.  Both  the  closure  constraint  and  the  P  current  con¬ 
straints  must  first  be  expressed  as  functions  of  7^(J^): 


M<l>q.  Zq) 


^  (2Zq) 


^^g;iii(^gsinc(27r7czg)  j^{k) 


[19a] 


J<p{^p  f  Zp) 


[19b] 


With  the  constraint  equations  written  as  above,  the 
functional  to  be  optimized  would  take  a  form  similar  to 


Eq.  [11]  and  the  optimal  j^{k)  obtained  as  before;  how¬ 
ever,  there  is  a  practical  issue  involved  in  the  constraint 
of  current  that  gives  reason  to  alter  Eq.  [19].  The  prescrip¬ 
tion  of  7^(0,  z)  at  a  finite  set  of  points  does  not  prevent 
large  values  of  /^((^,  z)  between  the  constraint  points.  As 
the  current  constraints  are  moved  closer  together  (requir¬ 
ing  a  larger  number  of  constraints  to  cover  the  same  area) 
non-zero  interconstraint  z)  values  correspond  to 
increasing  power  at  high  spatial  frequencies,  i.e.,  at  high 
k  values  in  j^[k].  Although  there  is  no  theoretical  limi¬ 
tation  on  the  spatial  frequency  content  of  /^((^>,  z),  any 
practical  design  will  involve  a  maximum  tolerable  spa¬ 
tial  frequency.  The  particular  frequency  will  depend  on 
the  application  and  will  often  be  determined  by  both  the 
width  and  minimum  bending  radius  of  the  conductor  to 
be  used  in  the  final  winding  of  the  coil.  For  these  reasons 
we  filter  out  the  impractically  high  frequency  compo¬ 
nents  of  /^((/),  z)  via  an  apodizing  function  w{k]: 


+»  ^ 

2  r”  [20] 

J77=— 00  —00 

where  in  the  present  work,  we  have  used  the  following 
functional  form  for  w(k): 


w(k)  =  1,  \k\<kt 


1 

2 


1  +  cos 


\k\  -  k\ 
kc-kj 


ic,<|ic|<Jc,  [21] 
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with  kf.  (“transition  frequency”)  and  (“cutoff  frequency”) 

chosen  based  on  the  desired  complexity  of  the  final  current 
density.  To  prevent  changes  in  the  direction  of  z)  from 
occurring  on  a  scale  smaller  than  Az  (for  example,  Az  could 
be  set  equal  to  twice  the  minimum  turning  radius  of  the 
wire  to  be  used  in  coil  construction),  k^  should  be  set  equal 
to  l/(2*Az);  kf.  is  then  set  to  allow  a  smooth  transition  in 
w[k]  while  preventing  the  attenuation  of  frequency  compo¬ 
nents  below  k^.  It  should  be  noted  that  the  spactug  of  the 
cmrent  constraints  is  set  equal  to  l/[2*kj  to  prevent  non¬ 
zero  values  of  current  density  between  constraiut  positions. 
Examples  of  specific  transition  and  cutoff  frequency  values 
used  in  practical  designs  will  be  given  in  the  Results  sec¬ 
tion.  We  now  apply  the  closure  and  current  constraints  to 
z)  given  above: 


+00  r 

■  ^  *■ 

m=-oo  ^  -00 


z,)=  2  I  clkr-^ik)j’S{k)  [22] 

r^ik)  =  i2Zg)e‘'"*«sinc{2iTkZg)w{k)  [23] 


and 


>  ^p) 


■il 


dk<t>^(k)j^ik)  [24] 


[25] 


The  apodizing  function  w(k]  must  also  be  incorporated 
into  the  constraint  expressions  to  assure  that  the  field 
constraints  will  be  met  by  the  constructed  current  den- 
sity: 


B 


+  00  ^  ^ 

>  ^n)  I 

m  =  — cc  —00 


dk^^ik)j^(k)  [26] 


iC;„(|27rica|)w(i:) 


[27] 


Only  the  three  sets  of  constraint  equations  are  modi¬ 
fied  with  the  inclusion  of  w{k);  the  inductance  expres¬ 
sion  of  Eq.  [12]  does  not  require  apodization.  The  inte¬ 
grand  of  Eq.  [12]  naturally  attenuates  high  spatial 
frequencies  more  extremely  than  w[k)  and  therefore  its 
inclusion  is  unnecessary.  The  attenuating  property  of  the 
inductance  integrand  can  be  seen  by  noting  that  for  large 
k,  I^{k)  ~  e^/(j27Tk],  K^J^k)  e~^/[j27rk),  and  therefore 
the  product  is  approximately  1/i:. 

With  the  apodization  discussed  above,  the  functional 
to  be  minimized  takes  the  form: 


U{i^m  =  Uj^rn  +  2  A„[B.(r„,  zj  - 


2  Ap[/^((^p,  Zp)  /(fip]  Ag[A(<]>q,  Zq)  Ag] 
p=l 


[28] 


where  the  inductance  expression  is  given  by  Eq.  [12]  and 
the  summation  over  all  constraints  is  that  of  Eq.  [11]  with 
N+P-l-1  constraints  based  on  Eqs.  [22],  [24],  and  [26]. 
Taking  the  derivative  of  U{j^[k]]  with  respect  to  /^(k), 
setting  the  result  equal  to  zero,  and  solving  fory^(k),  we 
get: 


mk)  =  2  +  E  Xpb^ik)  +  KqC^{k)  [29] 

n=l  p=l 

where 

a-^ik)  =  (— ]  w(k)  [30] 
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b^{k)  =  ( [31] 
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The  set  of  [A^,  A^,  A^}  (representing  AT+P+l  total  un¬ 
knowns)  is  found  by  substituting  j^{k)  expressed  by  Eq. 
[29]  into  the  N+P+1  constraint  Eqs.  [22],  [24],  and  [26]  to 
get  the  following: 
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[35] 


This  set  of  equations  can  be  written  as  a  single  matrix 
equation,  the  solution  of  which  is  a  vector  of  the  required 
Lagrange  multipliers: 


MA  =  c 


[36] 
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In  block  form  Eq.  [36]  is: 


where  the  six  distinct  blocks  of  the  matrix  are  all  known 
functions  of  the  constraint  coordinates  and  are  listed  in 
Table  1.  The  block  is  identical  to  that  used  in  Eq.  [16] 
but  derived  using  a^[k]  given  by  Eq.  [30]  instead  of  Eq. 
[15],  and  the  elements  are  functions  of  the  field  com 
straint  coordinates  only.  If  the  azimuthal  current  con¬ 
straints  are  included,  the  resulting  iV+Pby  N+P  matrix  is 
composed  of  M^,  and  M22;  here  M22  is  only  a 
function  of  the  current  constraint  coordinates  whereas 
Mi2  is  a  function  of  both  the  field  and  current  coordi¬ 
nates.  The  inclusion  of  the  closure  constraint  increases 
the  dimension  of  the  matrix  by  1  (now  of  rank  AT+P+l) 
and  M33  depends  only  on  the  closure  coordinates.  The 
M23  block  includes  both  closure  and  current  coordinates 
while  Mi3  has  both  closure  and  field  coordinates.  The 
Lagrange  multiplier  vector  would  be  {A^}  for  only  field 
constraints,  [A^,  A^)  for  field  and  current,  [A^,,  A^}  for  field 
and  closure,  and  (A^,  A^,  A^}  for  all  three  types  of  con¬ 
straints.  The  constant  vector  on  the  right  hand  side  of  Eq. 
[37]  is  partitioned  in  an  analogous  way. 

With  the  above  mathematical  forms  in  place,  the  de¬ 
sign  algorithm  is  simply  the  following.  First  a  set  of 
constraints  is  chosen  that  will  assure  that  the  gradient 
coil  will  fall  within  specifications.  The  field  constraints 
determine  the  type  of  gradient  coil  (transverse  versus 
longitudinal),  the  gradient  strength,  and  the  size  and 
absolute  placement  of  the  region  of  uniformity.  Current 
and  closure  constraints  are  chosen  such  that  the  physical 
coil  will  be  of  the  required  size  and  the  region  of  gradient 
uniformity  is  positioned  as  desired  relative  to  the  current 
density.  Once  all  constraints  are  set,  the  elements  of 
matrix  M  are  computed  via  the  known  forms  in  Table  1 
and  the  constraint  vector  is  assembled.  Equation  [37]  is 
solved  for  the  Lagrange  multipliers  through  Gaussian 
elimination  or  an  alternative  technique.  The  multipliers 
are  used  in  Eq.  [29]  to  yield  the  transform  of  which 
gives  the  current  density  of  minimum  inductance  satis¬ 
fying  all  of  the  specified  design  constraints. 

RESULTS  AND  DISCUSSION 

The  theory  in  the  previous  section  has  been  implemented 
in  the  MATLAB  software  language  (version  4.2,  The 
Mathworks  Inc.,  Natick,  MA)  and  used  to  create  a  variety 
of  gradient  coil  designs,  representative  examples  of 
which  are  presented  in  the  following  paragraphs.  The 
azimuthal  current  density,  /^((/),  z),  is  the  final  result  of 
the  above  analysis;  however,  the  discrete  wire  pattern 
approximating  the  current  density  is  obtained  via  stream 
function  methods  described  elsewhere  (7).  The  entire 
calculation,  from  specification  of  all  constraints  to  deter¬ 
mination  of  z)  and  the  approximation  with  discrete 
wire  paths,  typically  takes  less  than  three  minutes  on  a 
Sun  Microsystems  ULTRA  SPARC  1/140  workstation. 
All  analytic  forms  involved  in  the  evaluation  of  the  ma¬ 


trix  elements  of  Table  1  and  subsequent  construction  of 
7^(<^,  z)  are  well  known  (exponential  and  modified  Bessel 
functions)  and  there  are  readily  available  and  accurate 
routines  within  the  MATLAB  package  for  their  evalua¬ 
tion.  Once  the  final  wire  pattern  is  obtained,  the  pattern 
for  the  entire  coil  is  discretized  into  on  the  order  of 
10,000  to  20,000  elements,  allowing  for  numerical  veri¬ 
fication  of  the  achieved  magnetic  fields  and  inductance 
using  the  Biot-Savart  equation  and  Neumann’s  formula 
(22),  respectively.  Gradient  efficiency  (17,  defined  as  gra¬ 
dient  strength  per  unit  current)  can  be  gained  at  the  cost 
of  increasing  coil  inductance  simply  by  increasing  the 
number  of  wires  used  to  sample  the  current  density. 
Efficiency  varies  linearly  with  wire  number  while  the 
inductance  increases  as  the  square  of  the  wire  number. 
Inductance  also  varies  as  the  fifth  power  of  coil  radius 
(for  constant  efficiency  and  wire  pattern  (16));  therefore, 
for  a  given  design,  rj  •  is  a  constant,  independent 

of  both  radius  and  wire  number.  We  use  this  quantity  as 
a  figure-of-merit  or  performance  measure  in  comparing 
different  designs. 

The  most  familiar  gradient  coil  application  is  the  es¬ 
tablishment  of  a  uniform  gradient  field  within  the  central 
region  of  a  current  carrying  cylinder.  To  demonstrate  the 
effectiveness  of  the  combination  of  closure  and  current 
constraints  in  limiting  the  extent  of  the  final  current 
density,  the  four  different  combinations  of  constraints 
described  in  the  theory  section:  field,  field  and  closure, 
field  and  current,  and  field,  closure,  and  current,  were 
applied  to  the  problem  of  designing  an  extremely  short 
aspect  ratio  (1:1,  coil  length: coil  diameter)  central  region 
of  interest  (ROI)  coil.  Field  constraints  must  always  be 
included  because  the  minimum  inductance  current  den¬ 
sity  not  constrained  to  produce  any  field  is  simply  zero 
everywhere.  It  has  been  suggested  in  the  literature  (5) 
that  such  extremely  short  coil  designs  are  not  possible 
with  target  field  methods;  this  example  serves  to  demon¬ 
strate  that  using  the  methods  developed  in  this  paper, 
very  short  aspect  ratio  coils  can  be  designed  in  a  straight¬ 
forward  manner.  An  ROI  of  length  0.75*a  was  chosen 
over  which  four  field  constraints  were  placed:  =  a/4, 
4>n  =  0,  =  [-3,  -1, 1,  3}*(a/8)  with  where 

is  the  transverse  gradient  strength  in  T/m.  The  design 
corresponding  to  those  specifications  is  shown  in  Figs, 
la  and  lb.  The  coil  has  an  aspect  ratio  of  over  2.5:1  with 
the  edge  of  the  ROI  located  more  than  2a  away  from  the 
edge  of  the  coil.  To  limit  the  aspect  ratio  to  1:1  while 
maintaining  the  size  of  the  central  ROI,  a  closure  con¬ 
straint  was  placed  at  the  desired  position  (zj  and  a  series 
of  current  constraints  placed  over  the  adjacent  region; 
specifically,  the  closure  constraint  was  A  =  0  for  =  0, 
Zg  =  a  and  the  current  constraints  were  =  0  at  =  0, 
Zp  =  ±{(Zg  -h  dz^i);  for  i  =  0:55  and  dz  =  (3a/55)}.  The 
apodization  function,  w(k),  was  of  the  form  given  in  Eq. 
[21]  with  kc  =  l/[2^dz)  and  k^  =  0.6*k^,  set  to  prevent 
non-zero  /^(</),  z)  between  constraint  points.  This  consti¬ 
tutes  117  constraints  in  total  (AT  =  4,  P  =  2*56  =  112,  Q  = 
1)  and  the  results  are  shown  in  Figs.  Ic-lh.  The  effects  of 
the  closure  and  current  constraints  applied  separately  are 
shown  in  Figs.  Ic-ld  and  le-lf,  respectively.  Figures 
Ic-lh  demonstrate  that  the  combination  of  closure  and 
current  constraints  is  a  sufficient  condition  for  the  con- 
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FIG.  1.  Middle  ROI,  normalized  J^(z)s  with  corresponding  half¬ 
quadrant  wire  paths  for:  (a-b)  field  constraints  only  (short  heavy 
bars);  (c^d)  field  and  closure  (long  heavy  bar)  constraints;  (e-f)  field 
and  current  (long  bars);  (g-h)  field,  current,  and  closure  con¬ 
straints.  All  constraints  are  placed  symmetrically  about  z  =  0. 

straint  of  gradient  coil  length  and  they  must  be  used 
together  to  be  effective  in  that  regard.  The  calculated 
gradient  efficiencies  and  inductances  for  the  respective 
designs  are  shown  in  Table  2  for  a  =  0.325  m  (approxi¬ 
mate  radius  for  a  full  body  gradient  coil).  As  expected, 
the  merit  figure  and  the  aspect  ratio  decrease  as  con¬ 
straints  are  added;  however,  access  to  the  usable  region  is 
significantly  improved.  For  the  1:1  design,  the  region  of 
20%  gradient  deviation  extends  over  a  cylinder  of  diam- 

Table  2 

Gradient  Design  Example  Results  Summary 

Coil  Parameters 


Efficiency  Inductance  Merit 

_ [mT/m/A]  [fxH]  riO~T 

Central  ROI  designs'^ 

Field  constraints  0.090 

Field  +  closure  constraints  0.074 

Field  +  current  constraints  0.070 

Field  +  closure  +  current  0.055 

Edge  ROI  designs^ 

Field  constraints  0.061 

Field  +  closure  constraints  0.057 

Field  +  current  constraints  0.053 

Field  +  closure  +  current  0.032 

®  Merit  ~  {r}*  in  S.l.  units. 

Both  design  sets  are  for  a  =  0.325  m. 


eter  0.20  m  and  length  0.30  m  at  the  coil  center,  begin¬ 
ning  0.20  m  from  the  coil  edge,  whereas  for  the  2.5:1 
design,  the  same  region  extends  over  a  cylinder  of  diam¬ 
eter  0.30  m  and  length  0.40  m  at  the  coil  center,  begin¬ 
ning  0.75  m  from  the  coil  edge.  Figure  2a  shows  calcu¬ 
lated  contours  of  20%  gradient  deviation  on  the  xz  plane 
for  both  the  2.5:1  and  1:1  aspect  ratio  designs.  It  can  be 
seen  that  the  uniformity  of  the  gradient  field  is  not  com¬ 
promised  in  a  major  way  through  the  addition  of  the 
current  constraints.  The  wire  pattern  for  the  1:1  design  is 
only  slightly  more  complex  than  the  2.5:1  design,  the 
only  potential  concern  is  the  increased  current  density  at 
the  coil  edge.  This  increased  density  is  common  to  con¬ 
strained  length  gradient  coil  designs,  and  is  typically 
only  two  to  three  times  the  maximum  density  of  the  pure 
minimum  inductance  designs.  Increasing  the  aspect  ratio 
results  in  a  decreased  current  density  at  the  coil  edge  and 
the  two  can  be  traded  off  against  each  other  in  a  con¬ 
trolled  manner. 

As  the  above  example  suggests,  the  ability  to  restrict 
current  density  to  a  specified  region  of  the  coil  allows  for 
the  control  of  the  position  of  the  ROI  with  respect  to  the 
physical  end  of  the  coil;  an  important  extension  of  this 
idea  is  the  design  of  gradient  coils  with  a  usable  region  of 
uniformity  arbitrarily  close  to  the  end  of  the  coil.  To 
design  such  an  edge  gradient  coil,  a  set  of  field  con¬ 
straints  was  placed  over  the  edge  ROI  just  as  was  done  in 
the  design  described  above,  A  second  set  of  field  con¬ 
straints  was  placed  at  the  opposite  end  of  the  coil,  forcing 
the  first  moment  of  the  current  density  to  be  zero  by 
symmetry,  thereby  guaranteeing  a  torque  balanced  coil 


FIG.  2.  Contour  plots  of  20%  deviation  from  uniformity  over  xz 
planes  for  designs  with  and  without  current  constraints,  (a)  shows 
20%  contours  for  the  central  ROI  designs  with  (solid  line)  and 
without  (dotted  line)  current  constraints.  The  solid  vertical  line 
indicates  the  physical  edge  of  the  coil  for  the  1:1  aspect  ratio 
design,  (b)  shows  20%  contours  over  one  half  of  the  xz  plane  for 
the  edge  ROI  designs  with  (solid  line)  and  without  (dotted  line) 
current  constraints.  The  vertical  line  indicates  the  position  of  the 
physical  edge  of  the  coil  for  the  current  constrained  design. 
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FIG.  3.  Edge  ROI,  normalized  x}^(z)s  with  corresponding  half- 
quadrant  wire  paths  for:  (a-b)  field  constraints  only  (short  heavy 
bars);  (c-d)  field  and  closure  (long  heavy  bar)  constraints;  (e-f)  field 
and  current  (long  bars);  (g-h)  field,  current,  and  closure  con¬ 
straints.  All  constraints  are  placed  symmetrically  about  z  =  0.  Note 
relative  position  of  closure,  current,  and  field  constraints  yielding 
the  edge  gradient  coil. 


(20).  No  field  constraints  were  specified  over  the  central 
region  because  that  area  will  not  be  utilized  and  unnec¬ 
essary  constraints  only  serve  to  increase  the  inductance 
of  the  design.  The  end  of  the  coil  was  moved  adjacent  to 
the  ROI  by  applying  closure  at  the  outer  edges  of  the  field 
constraint  sets  and  a  set  of  current  constraints  applied 
over  the  region  exterior  to  the  closure  positions.  As  an 
example,  four  field  constraints  were  placed  over  a  region 
extending  from  z  =  a  to  z  =  2a:  =  a/4,  <j>^  =  0,  z^  =  {a  + 

(a/3)*i},  where  i  =  0:3  with  =  Gx*  ^n-  The  same  set 
was  placed  over  the  corresponding  region  on  the  negative 
side  of  the  coil,  constituting  eight  field  constraints  in 
total.  The  minimum  inductance  design  subject  to  these 
field  constraints  is  shown  in  Figs.  3a  and  3b.  The  aspect 
ratio  of  the  design  is  more  than  4:1  and  the  outer  edge  of 
the  ROI  is  more  than  2a  away  from  the  edge  of  the  coil. 
The  addition  of  the  closure  constraint  at  =  0,  z^  = 
max(zj  =  2a,  and  the  set  of  current  constraints:  J^p  =  0 
at  (j)p  =  0,  Zp  =  ±{(Zg  dzn);  for  i  =  0:45  and  dz  = 
(2a/45)},  constitutes  101  constraints  in  total  (iV  =  8,  P  = 
2*46  =  92,  Q  =  1)  and  as  in  the  central  ROI  design 
example,  w[k]  was  as  shown  in  Eq.  [21]  with  = 
l/(2*dz)  and  kf.  =  O.G^Jc^;  this  results  in  the  edge  gradient 
coil  design  shown  in  Figs.  3g  and  3h.  For  a  =  0.325  m 
(approximate  body  coil  dimension)  the  region  of  20% 
gradient  uniformity  begins  7  cm  outside  the  edge  of  the 


coil  and  ends  32  cm  inside.  For  comparison.  Figs.  3c-3f 
show  the  results  of  applying  the  closure  and  current 
constraints  separately.  The  calculated  gradient  efficien¬ 
cies  and  inductances  are  listed  in  Table  2.  Figure  2b 
shows  the  calculated  contours  of  20%  gradient  deviation 
on  the  xz  plane  for  designs  with  and  without  the  above 
current  constraints;  again,  it  can  be  seen  that  the  current 
constraints  do  not  significantly  affect  the  achieved  gra¬ 
dient  homogeneity.  The  edge  coil  wire  pattern  is  not 
significantly  more  complicated  than  the  unconstrained 
version  of  Fig.  3b  nor  is  the  maximum  current  density 
(located  at  the  end  of  the  coil)  extreme,  being  slightly 
more  than  twice  its  maximum  value  in  the  unconstrained 
example.  Just  as  in  the  1:1  central  ROI  design  example, 
the  positioning  of  the  ROI  at  the  coil  edge  is  accompanied 
by  a  decrease  in  the  merit  figure.  The  achieved  merit 
figure  depends  on  the  size  and  position  of  the  specified 
ROI  with  respect  to  the  allowed  extent  of  the  current 
density.  These  can  be  traded  off  against  each  other  for 
any  particular  application. 

The  ability  to  constrain  the  position  of  the  imaging 
region  with  respect  to  the  physical  gradient  coil  with 
target  field  methods  should  represent  an  important  im¬ 
provement  in  gradient  coil  design  in  general.  The  mini¬ 
mum  inductance  target  field  method  already  developed 
possesses  the  advantage  that  it  is  an  inverse  solution  to 
the  design  problem.  It  is  a  very  numerically  efficient 
technique  that  takes  advantage  of  the  FFT  algorithm  (23) 
and  its  use  allows  the  designer  to  gain  an  understanding 
of  the  connection  between  spatial  frequency  components 
of  the  magnetic  field  and  current  density.  We  believe  that 
one  of  the  major  limitations  of  target  field  methods  has 
been  the  lack  of  direct  control  over  the  nature  of  the  final 
form  of  the  current  density.  The  theory  described  in  this 
paper  removes  this  limitation.  There  are  several  impor¬ 
tant  applications  of  this  method.  The  constraint  of  total 
coil  length  for  a  middle  ROI  coil  will  allow  the  optimal 
design  of  gradients  for  use  with  short-bore  magnets.  It  is 
also  useful  to  limit  gradient  coil  lengths  simply  for  ease 
of  construction  and  potential  reduction  of  associated 
costs.  Head  sized  gradient  coil  design  will  especially 
benefit  from  this  method.  The  positioning  of  the  ROI 
directly  adjacent  to  the  coil  edge  will  allow,  for  the  first 
time,  imaging  of  the  neck  region  with  a  head-sized  gra¬ 
dient  coil.  It  should  also  allow  the  imaging  of  the  brain 
region  with  the  subject’s  eyes  positioned  at  or  outside  the 
coil  edge.  This  would  be  useful  in  functional  imaging 
studies  where  visual  communication  with  the  subject  is 
necessary.  We  have  designed  a  head-sized  gradient  coil 
capable  of  imaging  with  efficiencies  of  0.4  mT/m/A  over 
a  region  adjacent  to  the  coil  edge  and  having  an  induc¬ 
tance  of  800  juiH;  we  will  report  on  the  details  of  this 
design  and  experimental  results  of  imaging  with  this  coil 
in  a  forthcoming  paper.  The  ability  to  easily  design  gra¬ 
dients  of  this  nature  should  make  possible  new  imaging 
methods  that  were  previously  considered  impractical. 
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